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SUMMARY 


The  current  available  information  concerning  liquid- injection  thrub*-- 
vector  control  from  the  standpoint  of  both  system  performance  and  analysis  is 


summarized  in  this  report.  A  detailed  description  cf  the  complex  flow  field, 
created  by  the  injection  of  a  liquid  into  a  supersonic  nozzle  is  given  on  the 
basis  of  both  experimental  information,  including  pressure  measurements  and 
visual  evidence,  and  the  theory  of  compressible  fluid  flow.  Also  included 
is  a  brief  discussion  of  the  important  individual  processes,  such  as  pene¬ 
tration,  atomization,  vaporization,  and  chemical,  reaction,  which  are  integral 
parts  of  the  overall  flow  interaction.  Overall  system  performance  and  the 
effects  of  the  many  injection  parameters  are  then  considered  and  illustrated 
with  experimental  data  from  a  number  of  sources.  The  final  portion  of  the 
report  deals  with  analytical  methods  for  predicting  the  performance  of  liquid- 
injection  TVC  systems.  In  particular,  a  theoretical  model  based  on  a  detailed 
flow  field  solution pwhdch  hi!  “been  under ..developmencrat  Aerojet,  l  is  described 
in  detail £  Because  of  a  lack  of  funding,  this  model  was  not  fully  completed. 
Thus,  the  la^t  section  contains  the  specific  tasks  that  will  be  undertaken  to 
complete  the  mbdel  when  funding  again  becomes  available. 
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I. 


INTRODUCTION 


The  objectives  of  this  report  are: 

a.  to  promote  a  better  understanding  of  the  complex  flow  interaction 
associated  with  liquid-injection  thrust-vector  control  (LITVC); 

b.  to  provide  information  that  can  be  utilized  for  preliminary  LITVC 
system  design; 

c.  to  indicate  the  present  status  of  analytical  efforts  that  have  been 
directed  towards  the  development  of  a  theoretical  model  for  performance  prediction. 

Experimental  data  from  all  known  sources  to  date  have  been  compiled  with  the  first 
two  objectives  in  mind.  These  data,  especially  those  of  a  more  fundamental  nature 
such  as  measured  pressure  distributions  and  schlieren  photographs,  have  been  used  to 
supplement  compressible  flow  theory  in  the  discussions  of  the  various  facets  of  the 
liquid-injection  flow  field.  Actual  performance  data  have  been  carefully  selected  to 
illustrate  the  effect  of  each  injection  or  motor  parameter  on  system  performance.  This 
information,  as  presented  in  Section  III,  should  enable  the  design  engineer  to  deter¬ 
mine  approximately  the  performance  to  be  expected  from  any  given  system  and  should 
also  provide  a  means  for  selecting  the  proper  parameter  values  for  specific  applica¬ 
tions. 

The  sections  of  the  report  devoted  to  the  analytical  aspects  of  LITVC  have 
been  included  to  describe  the  approach  taken  here  at  Aerojet  to  develop  a  suitable 
comprehensive  model,  the  progress  that  has  been  made  to  date  on  that  model,  and  the 
steps  required  for  its  final  completion.  It  is  anticipated  that  soon  after  the  con¬ 
clusion  of  a  present  independent  research  and  development  experimental  program, 
which  will  supply  much  needed  flow  field  data,  another  report  will  follow;  this 
report  will  describe  the  finished  model. 
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I,  Introduction  (cont.) 
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The  feasibility  of  thrust-vector  control  by  fluid  injection  was  first  demon¬ 
strated  in  1952  by  G.  P.  Hausmann  at  United  Aircraft. ^  By  injecting  a  secondary  gas 
cross-stream  into  the  divergent  section  of  a  supersonic  nozzle,  Hausmann  showed  that 
a  lateral  force  exceeding  the  magnitude  of  the  secondary  stream  momentum  was  produced. 
The  additional  force  was  a  direct  result  of  the  interaction  of  the  two  streams  within 
the  nozzle,  which  caused  the  formation  of  a  local  high-pressure  region  along  the 

nozzle  wal3,;  in  the  vicinity  of  the  injector.  In  subsequent  tests  by  a  number  of 
2  3  and  H 

.Investigators ,  ’’  it  was  found  that  this  interaction  force  was  strongly  depend¬ 

ent  on  such  injection  parameters  as  the  location  and  configuration  of  the  injector, 
the  injectant  flow  rate,  and  the  properties  of  the  injectant  itself.  These  initial 
experiments  were  all  performed  with  gaseous  injectants,  but  because  of  material  prob¬ 
lems  with  hot  gases  and  bulky  system  requirements  for  cold  gases ,  the  use  of  a  liouid 
injectant  was  -proposed,  which  was  proved  feasible  by  NOTS  in  1959- 


'  The  first  class  of  liquid  injectants  to  be  considered  were  those  which  readily 
vaporized  when  encountering  the  high  temperatures  within  the  nozzle.  Of  all  the 
.liquids  of  this  type  tested,  bhe  Freon  group  emerged  as  the  most  successful.  However, 
it  Was  soon  discovered  that  injectants  which  chemically  react  with  the  main  stream 
gases  to  liberate  heat  would  produce  greater  side  forces  per  pound  of  injectant.  This 
-  led  to  the  investigation  of  such  liquids  as  nitrogen  tetroxide,  UDMH,  hydrazine, 
bromine,  and  perchloroethylene ,  which,  in  addition  to  giving  relatively  high  perform¬ 
ance,  also  created  toxicity  and  corrosion  problems.  To  combat  these  latter  effects, 

/  sodium  perchlorate  and  strcncium  perchlorate  solutions  in  water  were  tested  at 

Aerojet0’  and  were  found  to  have  completely  acceptable  storage  and  handling  charac¬ 
teristics,  while  still  producing  higher  performance  than  the  Freons. 

Initially,  the  injectors  consisted  only  of  a  single  port  in  the  nozzle  wall, 
but  further  efforts  to  increase  performance,  resulted  in  multiport  configurations  in 
which  the  ports  were  spread  circumferentially.  This  spreading  of  the  disturbance 
along  the  nozzle  wall  proved  to  be  successful  as  long  as  the  distance  between  adjacent 
ports  was  not  excessive. 
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1,  Introduction  (cont.) 

The  liquid-injection  TVC  concept  emerged  from  the  feasibility  stage  when  it 
was  adopted  for  use  on  three  major  solid-rocket  motors.  Both  the  Polaris  A3  and 
Minuteman  Wing  VI  second  stages  have  flown  successfully  with  an  injection  system 
using  Freon  114B2  and  a  three-port  Injector.  For  the  Titan  III  120-in. -dia  solid 
rocket  booster,  an  system  consisting  of  six  ports  evenly  spread  over  each  90° 

quadrant  was  selected.  Excellent  test  results  have  been  reported  for  this  configura¬ 
tion. 


The  liquid-injection  system  has  a  number  of  advantages  over  other  methods  of 
thrust-vector  control.  Since  there  are  no  mechanical  parts  exposed  to  the  propellant 
gases,  the"  material  problems  usually  associated  with  today's  high  temperature,  particle¬ 
laden  propellants  are  minimized.  Also,  the  injection  of  a  secondary  fluid  causes  an 
increase  in  the  axial  thrust  instead  of  a  thrust  loss  as  experienced  with  jetavators, 
jet  tabs,  and  swiveled  nozzles.  Another  very  useful  characteristic  of  this  type  of 
system  is  a  very  fast  response,  or  an  extremely  short  time  lag  between  the  command 
signal  and  application  of  the  full  control  force.  Since  there  is  no  break  in  the  flow 
passage  from  the  chamber  through  the  nozzle,  as  with  gimbaled  or  swiveled  nozzles, 
nozzle  design  is  simplified,  f  time  additional  advantages  include  applicability  to  all 
nozzle  sizes  and  a  decreasing  sy_ tern  weight  throughout  the  flight  as  the  injectanc.  is 
expended. 

Of  course,  this  system  also  has  certain  disadvantages.  T  im  a  performance 

/ 

standpoint,  the  most  serious  drawback  is  the  relatively  poor  performance  at  high 
jet-deflection  angles.  Deflections  over  6°  are  usually  difficult  to  obtain,  which 
limits  booster  applications.  The  entire  injection  system  car.  become  quite  complicated 
when  such  areas  as  pressurization,  tankage,  and  valving  are  carefully  considered. 
Finally,  tne  possibility  always  exists  of  the  injectant  supply  being  exhausted  before 
the  end  of  the  flight,  leaving  the  vehicle  without  means  of  control. 

The  suitability  of  this  system  for  a  particular  application  must  therefore  be 
del- rmined  through  trade-off  studies  considering  system  performance,  weight,  and 
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I,  Introduction  (cont.) 

reliability.  In  such  a  study  conducted  by  Space  General  for  a  typical  vehicle,  liquid 
injection  was  not  recommended  for  the  first-stage  motor,  but  was  found  to  be  very 
competitive  with  other  methods  of  thrust-vector  control  for  the  upper  stages. 5^- 
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II .  DISCUSSION  OF  FLOW  FIELD 

A.  GENERAL  DESCRIPTION 


When  a  liquid  is  injected  cross-stream  into  a  nozzle,  it  will  penetrate 
»  as  a  solid  jet  a  short  distance  into  the  supersonic  gas  stream  before  shattering  into 
droplets.  The  drag  forces  acting  upon  these  droplets  carry  them  downstream,  while- 
the  high  temper atures  of  the  gases  cause  the  injectant  to  vaporize.  The  vapor  thus 
generated  forms  a  body  along  the  nozzle  wall,  which  acts  as  an  obstruction  to  the 
primary  flow.  The  resulting  flow  field,  shown  schematically  in  Figure  1  for  a  single 
port,  is  very  similar  to  the  interference  pattern  produced  by  a  solid  blunt  body. 


The  primary  feature  of  this  flow  field  is  a  three-dimensional  bow  shock 
which  stands  slightly  upstream  of  the  injection  port  and  spreads  both  radially  and 
circumferentially  as  it  extends  downstream.  At  the  intersection  of  t;  e  shock  with  the 
nozzle  wall,  the  disturbance  can  spread  upstream  through  the  subsonic  portions  of  the 
boundary  layer.  If  the  shock  is  sufficiently  strong,  the  abrupt  streamwise  pressure 
gradient  can  cause  the  boundary  layer  to  separate  from  the  nozzle  wall. 


In  the  vicinity  of  the  injector,  where  the  shock  is  strongest,  the  boundary 
layer  separates  upstream  of  the  bow  shock  and  reattaches  to  the  vapor  body,  thus  modi¬ 
fying  the  effective  shape  of  the  obstacle.  The  presence  of  this  separation  region  at 
the  nose  of  the  vapor  body  causes  the  shocli  pattern  near  the  injector  to  differ  some¬ 
what  from  the  pattern  without  a  boundary  layer.  The  bow  shock  is  closer  to  the  body, 
and  the  compression  waves  generated  near  tne  separation  point  merge  into  an  oblique 
shock,  which  then  intersects  the  bow  shock  to  form  the  bifurcated  pattern  shown  in 
Figure  2a.  The  size  of  this  separation  region  depends  on  both  the  height  and  width 
of  the  body.  At  low  injectant  flow  rates  with  little  penetration,  the  separated  region 
merges  rather  smoothly  into  the  vapor  body,  causing  the  separation  shock  to  be  the  main 
shock  as  shown  in  Figure  3a.  At  the  other  extreme,  such  as  a  single  liquid  jet  pene¬ 
trating  far  into  the  nozzle  as  in  Figure  3b,  the  separation  region  will  be  relatively 
small,  with  the  bow  shock  the  dominating  feature  of  the  flow  field.  In  addition,  the 

size  of  the  separation  zone  for  a  given  body  increases  for  greater  boundary  layer 

•  > 

thicknesses. 
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II,  A,  General  Description  (eont.) 


Downstream  of  the  injector,  the  separated  boundary  layer  reattaches  to 
the  nozzle  wall,  as  shown  in  Figure  2b,  in  a  manner  similar  to  the  shock-boundary 
layer  interaction  produced  by  the  reflection  of  an  oblique  shock.  The  size  of  this 

-»  ~  *  r  „ 

separation  region  depends  on  the  shock  strength,  both  decreasing  as  the  shock  moves 


closer-  to  the  exit  plane. 


Photographic  evidence  of  some  of  these  flow  patterns  is  presented  in 
Figures  4  through  7-  The  shadowgraphs  shown  in  Figure  4  were  taken  at  the  JPL  20-in. 


Supersonic  Wind  Tunnel  during  the  injection  of  liquid  nitrogen  through  a  single  port 
*  _  .  _  l 3  t*. ^*UF.  ^  UnnV  /-.-P  O  m  9  nrtrl  Tint.h  vannr 


in  a  flat  -elate. ^  For  free  stream  Mach  numbers  of  2.01,  2.6l,  and  3.26,  both  vapor 
tody  and  -bow  shock  are  plainly  visible.  The  relatively  small  separation  region  can  be 
attributed  to  the  thin  boundary  layer  and.  most  probably ,  to  a  body  width  much  smaller 
than  the  body  height.  This  last  statement  can  only  be  postulated,  of  course,  since 
there  is  no  direct  evidence  of  the  vapor  "body  dimension  in  the  plane  of  the  plate. 
However,  it  does  appear  reasonable  because  of  the  relatively  small  diameter  of  the 
single  port,  the  high  injection  pressure,  and  small  .in jectant  flow  rate. 


Figure  5  is  a  schlieren  photograph  taken  at  the  same  JPL  wind  tunnel  for 

9 

M  =  2.01,  but  with  Freon  12  as  the  injectant  and  a  3-pcrt  injector-  In  this  case, 
the.  secaration  region  and  shock  are  more  prominent  than  in  the  .Liquid  nitrogen  case 
due  to  the  greater  lateral  spreading  of  the  vapor,  relatively  small  penetration,  and 


a  thicker  boundary  layer. 


The  effect  of  an  accelerating  free  stream  was  included  in  the  tests  per¬ 


formed  at  Stanford  with  injection  through  a  single  port  into  a  two-dimensional 
.  10  _ ..n  _ _ f...m/v  C.  P.m  1  -i fni'  fl  ni+vrtffott  1  c  nlli-hp  RT 


nozzle.10  The  overall  flow  pattern.  Figure  6,  for  liquid  nitrogen  is  quite  similar 
to  the  flat  plate  case.  In  addition  to  the  schlieren  photographs,  an  indication  of  the 
shock  and  body  dimensions  in  the  plane  of  the  nozzle  wall  was  also  found  In  these 
experiments  by  using  a  china-clay  technique  to  show  the  flow  pattern  along  the  wall. 
These  dimensions  proved  to  be  much  less  than  the  corresponding  ones  in  the  plane  of 
the  schlieren  photographs,  thus  supporting  the  presence  of  the  relatively  small  separa¬ 
tion  region. 
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II,  A,  General  Description  (cont.) 


For  an  axi symmetric  nozzle,  flow  visualization  techniques  must  he  replaced 
by  thorough  wall-pressure  measurements ,  to  determine  details  of  the  flow  field.  How¬ 
ever,  the  Minuteman  project  at  Aerojet,  during  an  actual  firing,  moused  a  camera  just 
outside  and  above  the  exit  plane  of  a  nozzle,  looking  down  the  inside  of  the  nozzle 
towards  the  ihro.rt.  Figure  7  is  one  of  the  photographs  taken  during  the  injection  of 
Freon  ll432  through  a  single  port,  and  shows  the  vapor  body  and  the  shock  along  the 
nozzle  wall. 


A  typical  pressure  distribution  for  Freon  llUBg  injection  through  a  3-port 
injector  in  the  Minuteman  Wing  VI  nozzle  is  shown  in  Figure  8.  This  nozzle  is  con¬ 
toured  with  an  overall  24.8:1  expansion  ratio,  arid  the  injector  is  located  at  an 
expansion  ratio  of  8.25:1.  This  pressure  distribution  is  typical  of  most  liquid- 
injection  cases  with  the  small  high-pressure  region  near  the  injector,  the  overall 
decrease  in  pressure  in  the  axial  direction  approaching  the  undisturbed  value,  and  the 
pressure  increase  from  the  injector  centerline  to  the  shock  in  the  circumferential 
direction  downstream  of  the  injector.  The  major  portion  of  the  total  side  force, 
about  70$,  is  contributed  by  the  downstream  region,  while  the  high  pressure  separation 
region  accounts  for  20$  and  the  injectant  momentum  the  final  10$. 

4o 

Similar  pressure  distributions  were  obtained  at  JPL  for  Freon  12  injected 
through  a  single  port  into  a  conical  nozzle.  With  the  single  port,  the  side-force 
contribution  from  the  separation  region  was  only  approximately  10$. 

B.  EFFECT  OF  VAPOR  BODY  SIZE  AHD  SHAPE 

The  formation  of  the  vapor  body  is  an  extremely  complex  phenomenon  in 
which  the  interrelated  processes  of  atomization,  vaporization,  mixing,  and,  in  seme 
cases,  reaction  occur  simultaneously.  The  size  and  shape  of  the  vapor  body  is  deter¬ 
mined  primarily  by  the  initial  penetration  of  the  liquid  stream,  the  rate  of  vapor 
generation,  and  the  injector  configuration.  To  produce  the  maximum  side  force,  the 
body  height  should  be  large  initially  to  produce  a  strong  shock  and  large  separation 
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IX ,  B,  Effect  of  Vapor  Body  Size  and  Shape  (cont.) 


X*egion,  and  continue  to  increase  downstream  to  prevent  the  pressure  levels  in  this 
relatively  large  area  from  falling  too  low.  A  stronger  shock  is  also  produced  by  a 
greater  body  width  as  the  three-dimensional  flow  over  and  around  the  body  approaches 
„  the  Uniting  two-dimensional  case  where  all  the  flow  must  go  over  the  body. 


4 


4, 


& 


The  initial  height  depends  upon  the  combination  of  penetration  and  rate 
of  vapor  generation.  For  low  penetration,  the  boundary  of  the  vapor  body  is  determined 
primarily  by  the  volume  of  vapor  generated,  which  is  then  free  to  spread  in  the  cross- 
sectional  plane  both  away  from  the  wall  and  along  it.  When  the  penetration  is  high, 
however,  the  body  height  is  mainly  a  function  of  the  trajectory  of  the  liquid  jet  and 
droplets,  with  a  corresponding  elongation  of  the  vapor  cross-section.  For  a  single 
port,  this  might  allow  more  of  the  primary  flow  to  spill  around  the  vfpor  body  rather 
than  being  forced  over  it,  which  would  then  result  in  lower  pressures  along  the  nozzle 
wall  away  from  the  body.  For  a  multiport  injector,  this  decrease  in  body  width  per 
port  might  permit  some  of  the  primary  flow  to  bleed  between  adjacent  ports,  also 
causing  lower  pressures.  This  initial  loss  could  be  compensated  for,  however,  if  the 
greater  mixing  provided  by  this  multiport  configuration  resulted  in  a  larger  vapor 
body  downstream  due  to  improved  vaporization  and/or  more  complete  chemical  reaction. 
Thus,  there  definitely  appears  to  be  an  optimum  penetration,  depending  upon  the 
in jectant  and  injector  configuration . 

After  this  initial  deflection  and  compression  of  the  primary  flow  caused  by 
the  "nose"  of  the  vapor  body,  the  flow  expands  again  as  it  follows  the  body  contour 
and  approaches  its  original  direction.  Thus  the  downstream  pressure  distribution  is 
governed  by  the  slope  of  the  vapor  body  boundary,  which  depends  on  the  amount  of 
vapor  generated,  the  mixing  of  the  vapor  and  free  stream,  and  the  degree  of  chemical 
reaction.  The  total  side  force  is,  therefore,  largely  determined  by  the  vaporization 
process. 


For  a  relatively  inert  in jectant,  such  as  Freon,  the  maximum  side  force  may 
be  produced  by  a  controlled,  rather  than  instantaneous  vaporization,  which  would  result 
in  an  increasing  vapor  body  size  approaching  the  nozzle  exit.  If  vaporization  occurs 
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II,  B,  Effect  of  Vapor  Body  Size  and  Shape  (cont.) 


too  quickly,  the  downstream  pressure  will  fall  rapidly  and  may  drop  below  the  undisturbed 
'■  j.iue  as  with  gas  injection.  On  the  other  hand,  the  maximum  body  size  will  r>t  be 
produced  if  the  vaporization  is  not  completed  before  reaching  the  nozzle  exit. 

For  a  reactive  injectant,  however,  vaporization  and  mixing  must  first  take 
place  before  the  chemical  reaction  can  be  initiated.  Thus  it  would  be  advantageous  in 
this  case  to  vaporize  the  injectant  quickly  to  give  the  vapor  a  chance  to  mix  and  react. 

C.  DETAILS  OF  INDIVIDUAL  PROCESSES 

1.  Penetration 


There  are  essentially  three  factors  affecting  the  penetration  of  a 
liquid  jet  into  a  cross-current  high  velocity  gas  stream: 

a.  the  momentum  of  the  jet  relative  to  the  free-stream; 

b.  the  break-up  of  the  jet  into  droplets; 

e.  the  shape  and  solidity  of  the  jet. 

Fcr  a  supersonic  free-stream,  additional  factors  may  enter  due  to 
the  shock-wave  formation  upstream  of  the  injector.  .Also,  vaporization  of  the  liquid 
and  possible  reaction  in  both  liquid  and  vapor  phares  could  produce  significant  effects. 

The  important  penetration  parameters  can  be  found  by  considering 
their  effect  on  the  above  factor-.  These  parameters  can  be  grouped  generally  as  free- 
stream  properties,  injectant  properties,  and  injector  design. 

The  primary  local  free-stream  parameters  are  static  and  dynamic 
pressures,  which,  for  a  given  gas,  are  functions  of  total  pressure  and  Mach  number. 

The  combination  of  low  total  pressure  and  high  Mach  number  will  produce  low  values  for 
coth  local  pressures,  thus  promoting  penetration.  The  secondary  effects  of  shocks  and 
vaporization  depend  on  the  free-stream  total  temperature  and  Mach  number,  while  the  gas 
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II,  C,  Details  of  individual  Processes  (cont.) 

composition  affects  any  chemical  reaction.  A  slight  reduction  in  penetration  would 
most  probably  result  from  these  secondary  effects. 

Finally,  the  injector  design  can  be  quite  significant  since  the  velocity, 
and  thus  momentum,  of  the  liquid  jet  is  determined  by  the  pressure  drop  across  the 
injector, and  the  port-discharge  coefficient.  The  penetration  is  also  a  direct  function 
of  the  port  size,  and  the  exact  form  of  this  relation  will  determine  whether  a  scaling 
effect  exists  for  any  particular  case.  Since  a  solid  stream  will  penetrate  farther 
than  a  spray,  the  jet  solidity,  which  depends  upon  the  contour  of  the  port  internal 
flow  passage,  is  also  an  important  factor. 

As  mentioned  previously,  there  does  appear  to  be  an  optimum  amount  of 
penetration  that  will  produce  the  greatest  side  force  for  any  particular  case.  In 
general,  the  optimum  value  should  be  greater  for  reactive  injectants  and  multi-port 
injectors.  If  a  scaling  effect  does  exist,  it  will  probably  be  more  difficult  to 
attain  this  best  penetration  in  the  larger  nozzles. 

2.  Atomization 


After  the  initial  penetration  of  the  liquid  jet,  the  stream  is  broken 

up  into  droplets  by  the  interaction  of  aerodynamic,  viscous,  and  surface  tension 

11  12  13  14 

forces.  Experimental  investigations  of  this  atomization  process  ’  ’  have  shewn 

that  the  size  of  these  droplets  depends  on  both  free  stream  and  injectant  properties, 
as  well  as  the  shape  of  the  jet.  To  promote  the  formation  of  smaller  droplets,  the 
free-stream  density,  and  the  velocity  should  be  high  and  the  injectant  viscosity  and 
surface  tension  low.  In  addition,  the  size  of  the  jet,  which  is  related  to  the  port 
diameter  and  injector  configuration,  should  be  small. 


Small  droplets  are  desirable  since  they  result  in  greater  amounts  of 
vapor  being  generated  within  the  nozzle.  This  can  be  attributed  mainly  to  the  addi¬ 
tional  total  surface  area  furnished  by  a  large  number  of  small  droplets  as  opposed  to 
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a  smaller  number  of  larger  droplets  comprising  the  same  total  mass.  This  larger  area 
•  allows  better  heat  transfer  to  the  liquid,  ar.d  thus  faster  vaporization.  '  ■: 


Two  different  types  of  vaporization  can  occur  within  the  nozzle: 
flash  vaporization  and  vaporization  due  to  heat  transfer.  Flash  vaporization  takes 
place  when  the  liquid  inject ant  is  suddenly  brought  into  an  environment  in  which  the 
pressure  is  below  the  vapor  pressure  of  the  liquid  at  its  initial  temperature.  The 
temperature  thus  must  drop  for  the  injectant  to  reach  an  equilibrium  state  at  the 
given  pressure  (see  Figure  9),  with  the  heat  loss  by  the  liquid  utilized  for  vapor 
generation.  The  amount  of  vapor  produced  in  this  manner  will  be  promoted  by  a  high 
initial  injectant  temperature,  a  low  final  pressure,  and  a  relatively  flat  vapor 
pressure  curve.  In  most  instances,  the  final  product  will  be  a  liquid- vapor  mixture, 
with  the  end  state  on  the  vapor  pressure  curve.  However,  in  extreme  cases,  the  liquid 
might  be  completely  vaporized,  so  that  the  final  state  will  be  in  the  vapor  region, 
and  thus  at  a  lower  temperature  than  the  liquid-vapor  mixture  at  the  same  pressure. 


The  most  prevalent  form  of  vaporization  ’is  a  direct  result  of  the 
heat  transfer  between  the  injectant  and  the  hot  gases  within  the  nozzle.  It  depends 
not  only  on  the  injectant  and  free-stream  properties  but  also  the  physical  mechanisms 
of  mixing  and  atomization.  In  general,  this  vaporization  will  be  promoted  by: 


a.  High  free-stream  temperatures; 

b.  High  free-stream  and  vapor  thermal  conductivity; 

c.  Low  free-stream  and  vapor  specific  heat; 

d.  Low  injectant  temperature; 

e.  Low  injectant  heat  of  vaporization; 

f.  Low  injectant  density; 

g.  High  relative  velocity  between  the  free  stream  and  liquid  droplets 

h.  Good  atomization  (small  droplets); 

1.  Good  mixing  between  the  free  stream  and  vapor. 
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XI,  C,  Details  of  Individual  Processes  (cont.) 

For  injection  into  sea-level  nozzles,  this  latter  form  of  vaporiza¬ 
tion  will  be  predominant,  but  flash  vaporization  could  become  significant  in  high- 
expansion-ratio  nozzles,  and  will  be  the  primary  form  in  cold-flow  experiments. 

U.  Reaction  and  Mixing 

If  an  exothermic  chemical  reaction  occurs  between  the  species  of 
the  free  stream  and  the  injectant  (liquid  or  vapor),  the  heat  generated  increases 
the  amount  of  jet  deflection  by : 

a.  causing  the  vapor  body  to  expand,  thus  forming  a  larger 
obstruction i 

<'* 

b.  transmitting  heat  into  the  primary  flow,  which  has  been 
shown-*- 5  to  be  analogous  to  mass  addition  with  resulting 
higher  pressures. 

The  molecular  weight  of  the  reaction  products  as  compared  to  that  of  the  reactants  is 
also  important  in  determining  the  final  body  size.  The  choice  of  an  optimum  reactive 
injectant  will  thus  depend  on  the  composition  of  the  primary  flow,  the  heat  of  reaction, 
and  the  reaction  products.  To  encourage  this  reaction, ~the  mixing  of  the  injectant  with  _ 
the  free  stream  should  be  promoted,  and  since  most  reactions  occur  with  the  injectant 
vapor  phase,  vaporization  should  be  accelerated. 

As  more  of  the  free-stream  gases  mix  with  the  injectant  vapor,  two 
opposing  effects  occur.  The  vapor -body  size  is  increased  due  to  higher  temperatures, 
greater  vaporization,  and  greater  reaction,  but  less  of  the  free  stream  is  being 
deflected.  However,  in  a  study  conducted  at  Aerojet  using  a  computer  program  for  tur- 
bulent  mixing  with  chemical  reactions,  °  the  overall  effect  of  greater  mixing  was  the 
production  of  a  larger  effective  vapor  body,  and  thus  greater  jet  deflection  or  side 
fores.  This  was  especially  true  in. the  case  of  chemical  reaction. 
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III.  PERFORMANCE  CHARACTERISTICS  v 

A.  TOPICAL  SYSTEM  PERFORMANCE  '' 

'  -  •»  '  O 

,  ^  - 

The  final  criterion  to  be  used  in  evaluating  any  given  liquid-injection 

thrust-vector-control  system  is,  of  course,  the  performance  of  that  system,  i.e.  the 
amount  of  side  force  produced  by  a  given  injectant  flow  rate,  or  conversely,  the 
mjeetant  flow  rate  required  for  a  given  amount  of  side  force.  To  present  system 
performance,  either  of  the  curves  shown  in  Figure  10  are  most  frequently  used. 

The  first,  sometimes  referred  to_  as  the  ’’gain  curve",  shows  side  force  ~ 
as  a  function  of  injectant  flow  rate,  with  the  motor  thrust  and  mass  flow  used  to 
nondxmensionalize  the  variables.  It  is  convenient  to  present  the  performance  in  terms 
of  these  ratios  since  the  side  force  ratio  (Fc/F^)  is  the  arctangent  of  the  resultant 
jet  deflection  angle,  and  the  mass  flow  ratio  is  a  scaling  parameter  that 

enables  the  performance  of  different  size  motors  to  be  compared. 

The  second  curve  shows  the  injectant  specific  impulse,  (ls)s,  or  -he  side 
force  per  unit  injectant  flow  rate,  as  a  function  of  the  side-force  ratio.  Since  the 
system  requirements  are  usually  stated  in  terms  of  jet  deflection,  this  latter  curve 
allows  the  efficiency  of  the  LlfVC  system  to  be  evaluated*  for  given  requirements. 


Either  of  these  two  curves  can  readily  be  transformed  to  the  other,  since. 

for  a  constant  motor  specific  impulse,  straight  lines  emanating  from  the  origin 

represent  lines  of  constant  (I  )s  in  the  gain  curve,  and  constant  w  /w.  in  the  I 

s  s  A  s 

curve . 

The  shape  of  the  curves  in  Figure  10  are  typical  for  a  constant  injection 
pressure  system.  As  the  injectant  flow  rate  increases,  the  slope  of  the  gain  curve 
tends  to  level  off,  and  m  some  cases,  will  actually  become  negative  due  to  the 
reflection  of  tne  shock  off  the  opposite  nozzle  wall.  Thus,  for  a  given  set  of 
injection  conditions,  there  is  a  maximum  side  force,  which  usually  occurs  at  a  very 
high  injectant  flow  rate. 
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III,  A,  Typical  System  Performance  (cont.) 


v  "  _  The- loss  in  system  efficiency  at  the  higher  mass  flow  and  side  force 

ratios  is  indicated  by  the  I  curve  of  Figure  10.  The  higher  the  jet-deflection 
requirement,  the  lower  the  injectant  specific  impulse.  Thus,  at  the  maximum  side- 
force  condition  for  the  system,  the  (I  )s  is  usually  so  low,  that  it  would  not  he 

v*;,  S 

desirable  to  operate  the  system  at  this  point. 


1  I 


B.  SIGNIFICANT  PARAMETERS 

■  •  Since  the  mechanism  producing  the  side  force  is  so  complex,  it  is  not 
surprising  to  find  that  system  performance  is  sensitive  to  a  larg^  number  of  para¬ 
meters.  These  parameters  can  be  grouped  in  the  following  manner: 

1.  Injection  parameters 

a.  injectant  properties 

b.  port  diameter  and  injection  pressure 

.  c.  injector  location 

d.  injector  configuration 

e.  angle  of  injection 


2.  Motor  parameters 

a.  chamber  conditions 

b.  propellant  properties 

'  “  c.  nozzle  size  and  shape 

In  general,  the  effect  of  each  parameter  can  be  deduced  from  experimental 

* 

data,  and  by  considering  the  flow  field  and  individual  processes  discussed  in  the 
previous  section,  these  effects  can  usually  be  logically  explained.  However,  with  so 
many  parameters,  it  is  sometimes  very  difficult  to  isolate  a  single  parameter  and  its 
effects,  since  a  change  in  one  parameter  may  necessarily  cause  a  corresponding  change 
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>r  „  I 

•  - 

in  some  other  parameter.  An  example  of  this  is  the  relationship  "between  injectant  flow 
rate,  port  area,  and  injection  pressure*  It  is  not  possible  to  change  only  one  of  these 
parameters  for  a  given  injectant.  - 

1.  Injectant  Properties  '  ‘ 


Liquid  injectancs  can  in  general  be  grouped  into  three  classes: 

a.  inert  liquids,  - 

b.  reactive  liquids, 

c.  monopropellants  or  bipropellants. 


The  inert  liquids  essentially  do  not  react  with  the  propellant  gases 
of  the  primary  flow,  and  thus  depend  primarily  upon  vaporization  to  produce  the  body 
of  vapor  which  forms  the  flow  obstruction.  Since  the  vaporization  requires  energy  from 
the  free  stream,  the  effective  flu-  deflection  is  somewhat  reduced.  In  the  inert 
liquid  class,  the  Freon  compounds  have  been  found  to  produce  the  best  performance  due 
largely  to  their  low  heats  of  vaporization.  Freon  113  and  llUBg  have  been  extensively 
used  since  their  relatively  high  liquid  density  reduces  tankage  requirements. 


Reactive  liquids,  due  to  chemical  reaction  with  the  propellant  gases, 
usually  liberate  heat  within  the  nozzle,  which  not  only  tends  to  increase  the  size  of 
the  vapor  body  by  increasing  temperature  and  the  amount  of  gas  generation,  but,  in 
transferring  energy  to  the  primary  flow,  also  tends  to  increase  the  flow  deflection. 

Since  most  propellant  gases  are  fuel  rich,  the  most  successful  reactive  injectants  have 
been  oxidizers,  such  as  nitrogen  tetroxide  and  hydrogen  peroxide.  Because  of  storability 
and  corrosion  problems  associated  with  these  particular  liquids,  a  new  class  of  liquid 
oxidizers  arose  -  perchlorate  salt  solutions.  Excellent  performance  has  been  obtained 
with  saturated  water  solutions  of  sodium,  strontium,  and  lead-perchlorates,  and  by 
using  a  reactive  solvent  such  as  hydrogen  peroxide,  this  performance  was  further 
increased.  With  a  reactive  solvent,  however,  the  storabij.ity  and  corrosion  problems 
again  arise,  but  would  not  be  as  significant  as  in  the  case  of  the  solvent  alone  used 
as  the  injectant. 
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In  the  third  class  of  Injectants,  the  monopropellants  and  bipropellants 
the  injectant  itself  reacts  as  soon  as  it  is  raised  to  a  sufficiently  high  temperature. 
Thus,  the  performance  of  these  injectants  does  not  greatly  depend  on  the  composition  cf 
the  propellant  gases  as  with  the  reactive  liquids.  In  actual  tests  with  typical  liquid 
and  solid  propellants,  the  performance  of  such  inject, ants  as  hydrazine,  UDMH,  Cavea-B 
nitromethane ,  UDMH-hydrazine ,  and  UDMH-IRFNA  have  in  general  been  greater  than  the 
inert  liquids,  but  less  than  the  oxidizers .  One  scheme  that  has  been  proposed  to 
increase  the  performance  of  monopropellants  and  bipropellants  is  the  use  of  a  pre¬ 
combustion  chamber  to  initiate  the  chemical  reaction  just  before  the  injectant  enters 
the  nozzle.  Thus,  some  of  the  greater  performance  associated  with  a  gas  injectant  could 
be  realized. 

hi  h2 

A  wide  range  of  liquid  injectants  have  been  tested  at  Aerojet  ’  , 

5  30 

HOTS  ,  and  Lockheed  in  actual  hot  firings  and  some  of  the  results  of  these  tests  are 

k 

shown  in  Figures  11C*  and  12.  In  these  figures ,  the  relative  performance  of  different 
injectants  can  be  compared  for  the  ame  injection  and  motor  parameters.  It  should  be 
noted,  however,  that  some  injectants  may  be  more  sensitive  than  others  to  certair 
parameters,  which  could  cause  a  shift  in  relative  performance  levels.  A  good  example 
of  this  is  the  greater  effect  the  propellant  gas  composition  and  temperature  have  on 
the  performance  of  reactive  liquids  than  on  inert  liquids,  monopropellants  and  bipropel¬ 
lants  .  , 

In  addition  to  performance,  the  final  selection  of  an  injectant  for 
a  given, application  will  depend  greatly  on  such  properties  as  density,  corrosiveness, 
toxicity,  and  storability,  which  will  affect  the  size,  weight,  complexity,  and  relia¬ 
bility  of  the  injection  system. 


*  Figure  numbers  designated  with  the  letter  "C"  are  contained  in  the  confidential 
supplement. to  this  report'.--.;  „  - 
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2.  Port  Diameter  and  Injection  Pressure 

The  effects  of  port  diameter  and  injection  pressure  cannot  he 

determined  separately,  since  at  a  given  in„ectant  flow  rate  one  cannot  he  varied 

without  causing  a  corresponding  change  in  the  other.  The  port  diameter  is  an 

important  parameter  in  the  atomization  process,  with  smaller  droplet  sizes,  and 

subsequently  faster  and  more  efficient  vaporization,  resulting  from  the  smaller 

port  sizes.  The  penetration  of  the  liquid  jet,  which  affects  hoth  the  initial 

strength-, of  the  main  shock  and  the  size  of  the  separation  region,  is  primarily  a 

function  of  the  initial  injectant  velocity,  which  in  turn'  depends  greatly  on  the 

injection  pressure.  Therefore,  in  most  cases  it  would  appear  that  the  comhination 

of  a  small  port  and  high  injection  pressure  would  produce  the  hest  performance 

at  a  given  injectant  flow  rate.  This  trend  is  illustrated  in  figure  13C,  where 

the  lines  of  constant  injection  pressure  and  constant  port  diameter  are  shown  in 

the  I  plot  for  a  single  port  injector, 
s 

The  question  that  now  naturally  arises  is  which  of  the  two  methods 
of  varying  injectant  flow  rate  is  preferable  -  constant  port  area  with  variable 
injection  pressure,  or  constant  injection  pressure  with  variable  port  area.  To 
compare  the  performance  of  the  two  systems,  they  must  he  matched  at  some  particular 
design  point.  The  natural  choice  for  this  design  point  is  the  maximum  injectant 
flew  rate,  for  which  the  highest  injection  pressure  and  largest  port  size  would 
be  required.  At  the  lower  flow  rates  then,  the  constant  pressure  system  would 
have  a  larger  injection  pressure  and  smaller  port  size,  thus  producing  greater 
penetration  and  higher  aide  forces,  lb.-,  performance  of  the  constant  area  system 
would  be  especially  poor  at  the  very  low  flow  rates,  where  with  the  low  injection  . 
pressures,  the  injectant  may  barely  penetrate  the  boundary  layer., 

A  good  comparison  of  the  two  systems  is  given  in  Figure  14C  where 
the  typical  performance  curves  for  each  system  are  shown.  Unlike  the  curves 
previously  discussed  for  a  constant-  pressure  system,  the  gain  curve  for  a  constant 


Page  17 


tm-i6-sro 


i  :  :«s  -  * 


,  ■'  j*  -*  < , 
C  f  •*-  “  *  * 
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area  system  has  an  inflection  point,  with  a  corresponding  maximum  in  the  X  curve.  V? 

>  '  »  8 

thus,  in  the  low  flow  range,  the  efficiency  of  a  constant  area  system  decreases 

at  the  smaller  injectant  flow  rates,  while  the  constant  pressure  system  efficiency  ’ 

increases . 

.  ^ 


She  trend  of  increasing  performance  with  increasing  injection 
pressure  at  a  given  flow  rate  does  not  held  indefinitely  as  illustrated  in  Figure 
15C.  In  the  previous  section,  it  was  concluded  that  an  optimum  amount  of  penetra¬ 
tion  did  exist.  Phis  conclusion  is  supported  by  the  fact  that  there  is  an  optimum 
injection  pressure,  which  varies  with  motor  chamber  pressure.  For  higher  pressures 
within  the  nozzle,  higher  injection  pressures  are  required  to  reach  the  optimum 
penetration . 

3.  Injector  Location 


There  are  a  great  number  of  factors  to  be  considered  in  determin¬ 
ing  the -optimum  injector  location..  From  an  aerodynamic  viewpoint,  as  the  injector 
is  moved  closer  to  the  throat,  the  shock  wave  becomes  steeper,  for  a  given  flow 

deflection,  because  of  the  lower  free-stream  Mach,  number.  Thus  the  nozzle  wall 

% 

area  covered  by  the  disturbance  increases.  Also,  despite  the  fact  that  the  shock 
strength  is  less  at  the  lower  Mach  numbers,  the  difference  between  the  pressures 
behind  the  shock  and  the  free-stream  pressures  is  still  greater  for  locations 
closer  to  the  throat  since  the  free-stream  pressures  are  higher.  The  combination 
of  larger  area  and  greater  pressure  differences  would  produce  higher  side  forces, 
so  that  from  strictly  aerodynamic  considerations,  the  optimum  injector  location 
should  be  as  close  to  the  throat  as  possible .  However,  two  other  important  factors 
must  be  considered  -  nozzle  geometry  and  vapor  body  size. 


Due  to  the  circular  cross-section  of  a  nozzle,  the  nozzle  geometry 
tends  to  reduce  performance  as  the  injector  is  moved  closer  to  the  throat.  This 


* 
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effect  can  best  be  explained  by  first  setting  up  a  cylindrical  coordinate  system 
*  (r,  9,  x)  with  the  nozzle  center  line  as  the  x  axis  and  the  injector  located  on  the 
nozzle  wall  at  0  =  0°. 


1 


Since  the  pressure  force  always  acts  normal  to  the  nozzle  wall,  a  differential  force 

dF  =  PdA  =  Prd  8-^— 
cos  a 

acting  at  any  point  on  the  wall,  will  have  the  following  component  contributing 
directly  to  the  side  force: 


dF  =  dF  cos  6  cos  a 
s 

=  (P  r  d9  dx)cos  0 

As  ©increases,  dfg  becomes  smaller,  and  for  90°<8<270°,  a  -- 

negative  side  force  component  is  produced.  Therefore,  as  the  pressure  disturbance 
spreads  along  the  nozzle  wall  in  the  c  rcumferential  direction,  it  becomes  less 
effective  in  producing  a  side  force.  «s  the  disturbance  spreads  substantially  into 
the  upper  half  of  the  nozzle,  an  actual  decrease  in  side  force  may  result. 

This  performance  loss  due  to  geometry,  sometimes  called  "curvature 
loss,"  will  become  more  predominant  as  the  injector  is  moved  closer  to  the  throat. 
The  steeper  shock  waves  at  these  locations  tend, to  further  aggravate  the  problem. 
Figure  16  illustrates  a  situation  in  which  the  injector  was  too  close  to  the  throat 
m  a  Mmuteman  nozzle  ^  and  the  shock  wave  was  hitting  the  opposite  nozzle  wall. 
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Pressure  taps  were  located  at  various  x  stations  at  Q  =  l80° ,  and  the  shock 
reflection  was  clearly  indicated  by  a  sharp  rise  in  pressure,  which  moved  upstream  s 
as  the  injeetant  flow  rate  increased.  The  gain  curve  for  this  test  is  shown  in 
Figure  2k ,  and  is  extremely  "flat”  for  w  /w,  >0.15. 

'  ^  S  ' 

* 

The  final  factor  that  must  be  considered  in  a  discussion  of  the 
optimum  injector  location  is  the  size  of  the  vapor  body,  which  determine.s  the  initial 
flow  deflection  and  the  downstream  pressure  distribution.  For  a  given  flow  rate  and 
injection  pressure,  the  penetration  of  the  liquid  jet  will  be  less  at  locations 
nearer  the  throat  and  the  volume  occupied  by  the  resulting  vapor  near  the  injector 
will  be  smaller  due  to  the  higher  free  stream  pressures.  Thus,  the  initial  flow 
deflection  will  be  less .  However,  the  higher  f ree-stream  temperatures  and  greater 
distance  to  the  nozzle  exit  will  promote  vaporization,  mixing,  and  reaction,  which 
will  tend  to  produce  higher  downstream  pressures  as  long  as  these  processes  have 
not  been  completed  within  the  nozzle . 

To  summarize  then,  the  injector  should  be  located 

* 

a.  Close  to  the  throat  to  optimize  the  aerodynamic  effects; 

b.  far  from  the  throat  to  minimize  the  losses  due  bo  the 
geometrical  effects; 

c.  at  some  intermediate  point  to  optimize  vapor  body  size. 

Thus,  there  is  definitely  some  optimum  location,  which  will  depend  upon  the 
relative  importance  of  each  of  the  above  three  factors .  Such  an  optimum  location 
is  clearly  shown  in  Figure  IT,  where  for  this  particular  case,  the  maximum  side 
force  for  a  given  flow  rate  would  be  produced  by  injection  at  x/L  =  0  J+5 
(x  =  axial  distance  from  throat  to  injector,  L  =  nozzle  length).'’  As  the  flow  rate 
is  increased,  the  location  of  this  maximum  usually  shifts  slightly  down  stream 
as  the  geometrical  effects  become  more  significant. 
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A  multiport  injector  also  tends  to  shift  the  optimum  location 
downstream,  since  this  injector  spreads  the  disturbance  more  circumferentially 
than  a  single  port.  The  curvature  losses,  therefore,  become  effective  sooner 
as  the  injector  is  moved  towards  the  throat.  Figure  18C  compares  the  effect 
of  injector  location  on  a  single  port  and  multiport  injector  for: the  same  nozzle t 

One  point  that  has  been  a  source  of  confusion  for  some  time  is 
the  variation  of  optimum  location  from  nozzle  to  nozzle .  This  variation  has  been 
quite  significant,  ranging  from  x/L  =  0.2  to  x/L  =0.6.  The  cause  of  such  an 
inconsistency  may  lie  once  again  with  the  geometrical  effects .  For  a  given  x/L, 
snock  reflection  and  curvature  losses  will  be,  more  of  a  problem  with  a  long ,  - , 

*  ",  _  v 

o lender  nozzle  than  with  a  short  stubby  one.  Thus,  the  ratio  of  nozzle  length  to 
exit  diameter,,  L/d,  appears  to  be  a  possible  correlating  parameter.  Figure  19C 
shows  that  such  a  correlation  does  exist,  with  the  .optimum  location  moving  down¬ 
stream  as  L/d  increases.  Of  course,  the  optimum  x/L  will  vary  somewhat  for  a 
given  nozzle,  depending  on  the  various  injection  parameters.,,  and  thus  only  the 
general  trend  can  be  indicated.  Also,  in  some  cases,  an  actual  optimum  was  not 
reached,  but  was  either  upstream  or  downstream  of  the  range  tested.  This  is 
indicated  by  an  arrow  in  Figure  19C.  . 


Experimental  data  indicate  that  increased  performance  can  be 

gained  by  injecting  the  liquid  upstream  against  the  flow  of  the  main  stream,  as 

can  be  seen  in  Figure  20.  One  of  the  main  effects  of  this  upstream  injection  * 

angle  is  an  increase  in" the  size  of  the  separation  region  upstream  of  the  injector. 

40 

Actual  pressure  distributions  measured  at  JPL  show  that  separation  was  initiated 

♦ 

more  than  twice  the  distance .upstream  of  the  injector  when  the  injection  angle  was 

t 

increased  from  0°  to  40°  (measured  with!  respect  to  the  .plane  normal  to  the  nozzle 
center  line,  with  upstream  values  positive) .  The  pressures  'downstream  of  the 
injector  were  found  to  be  slightly  higher  and  spread  more  circumferentially  with 
the  40°  angle.  -  •"  -  -  _  -  *,  •'  4  ? 
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The  larger  separation  region  with  an  upstream  injection  angle 
■  can  b’e  attributed  to  a  greater  penetration  with  respect  to  the  nozzle  wall,  as 
well  as  a  'slight  upstream  shift  in  the  actual  physical  obstruction'.  The  higher 
'relative  velocities. between  injectant  and  primary  flow  near  the  injector,  and 
longer  stay  time  for  the  injectant  within  the  nozzle,  promote  atomization, 
vaporization,  mixing,  and  reaction,  which  increase  the  vapor  body  size  and  thus 
account  for,  the  larger  downstream  pressure  area .  ,  .•  * 

*  For  cases  where  the  shock- reflection  effects  predominate,  (i.e. , 

high  flow -rate.s  and  locations  near  the  throat),  the  curvature  losses  with  the 

larger  upstrepi  angles  wiil  be  greater,  thus  reversing  the  usual  trend.  This 

situation^is  illustrat.ed  in  Figure  21,  which  shows^a  ,0°  injection  angle  to  be 

superior. to  a, 30°  angle  for  x/L  =  0.2,  but  the  30°  angle  to  be  better  farther  • 

3  3  '  34 

downstream..  3  Also,  the  effect  of  flow  rate  is.  evident 'since  the  crossing 

* 

of  the  0°  and  30°  performance  curves  occurs  farther  downstream  at  the  higher 

*  -  V  J 

flow  rates.'  *  .  ,  . 

t  »  -  ’ 

*  **  k  .  . 

5*  Injector  Configuration  -•  *  ", 

The  multiport  injector  usually  produces  higher  side  specific 
impulses  than  a  single  port  due  to  the  greater  spread  of  the  disturbance  area. 
Also,  by  distributing  the  injectant,  the  intermediate  processes  of  atomization, 
vaporization,  mixing  and  reaction  are  greatly  enhanced,  so  that  a  multiport 
injector  would  definitely  be  recommended  for  a  reactive  liquid.  The  only  port 
arrangement  -that  has  proved  successful  is  a  circumferential  spacing  at  the  same 
axial  coordinate .  Therefore,  it  is  necessary  to  consider  two  additional  variables 
with  multi-port  injectors-  -the  spacing  between  ports  and  the  number  of  ports . 

-  \  • 

Figure  22C  shows  that  there  is  an  optimum  spacing  for  a  given 

port  size  and  number  pf  ports.  When  the  ports  are  too  close  together,  mutual 
interference  occurs  between  the  disturbances  generated  from  each  port,  thereby. 
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reducing  the , effectiveness  of  the  interaction  process.  If  the  ports  are  too, far 
apart,  some  of  the  primary  flow  may  go  between  the  ports  essentially  unaffected, 
and  the  greater  spread  of  the  disturbance  will  once  again  bring  the  curvature  "  .. 
losses  into  effect.  •  - 

Available  data  indicate  that  increasing  the  number  of  ports  has 
a  much  greater  effect  on  reactive  liquids  than  on  inert  injectants .  This  can  be  •>, 
attributed  to  the  better  mixing  and  mo^e  complete  chemical  reaction  that  results 
from  the  greater  dispersion  of  the  injectant.  Figure  23C  shows  the 'increase  ih" 
tne  performance  of  a  reactive  liquid  as  the  number  of  ports  is  increased. ,  However, 
since  the  curvature  losses  are  also  increasing  as  the  total  spread  of  the  ports 
increases,  a  point  is  finally  reached  where  performance  will  drop  if  additional 
ports  are  added. 

The  best  combination  of  port  spacing  and  number  of  ports  will 
depend  strongly  on  the  particular  injectant,  and  will  also  vary  with  certain 
injection  parameters,  particularly  injector  location. 


In  discussing  the  injector,  one  factor,  that  should  be  mentioned 

is  the  shape  and  solidity  of  the  liquid  jet  as  it  enters  the  nozzle.  This  factor 

will  greatly  influence  the  penetration  of  the  jet,  as  well  as  the  atomization  process 

43 

to  a  lesser  extent.  Lockheed  stated  that  for  both  a  constant  area  orfice  and 
variable  area  valve,  the  liquid  stream  may  separate  from  the  flow  passage  walls, 
changing  a  solid  stream  to  a  spray like  jet,  as  shown  below. 


Constant  Area  Orifice 


Variable  Area  Value 


Without 

Separation 


With 

Separation 
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This  separation  is  promc*  1  by  a  short  flow  passage,  a  high  flow  rate,  or  a  low 
back  pressure,  and  will  greatly  depend  on  the  shape  of  the  passage.  In  tests 
with  a,  constant  area  orifice  at  Lockheed,  as  the  injection  pressure  was  increased, 
a  discontinuity  developed  in  the  flow  rate  because  of  the  onset  of  separation. 
Such  a  phenomenon  was  also  noticed  at  UTC  for  a  variable  area  valve,  and  could 
explain1  some  of  the  irregularity  of  liquid  injection  data. 


.  ,  v:h  < 

'21 
•  ?  -1 


This  separation  of  the  liquid  stream  causes  a  decrease  in  side 
force  since  the  jet  penetration  is  drastically  reduced.  This  effect  appears  to 
dominate  over  the  Detter  atomization,  which  would  improve  performance.  A  good 
example  of  the  separation  effect  was  the  loss  in  performance  that  occurred  with 
the  Minuteman  LITVC  system  when  the  shape  of  the  pintle  in  the  variable-area  val /e 
was  changed  to  give  better  flow  control.  Bench  tests  showed  that  the  new  valve 
produced  a  definite  spray  pattern,  in  contrast  to  a  solid  stream  from  the  old  valve . 

It  thus  appears  that  the  performance  of  any  LITVC  system  can  be 
directly  affected  by  the  flow  characteristics "of  the  liquid  stream,  and  the  evi¬ 
dence  to  date  indicates  a  solid  rod-like  flow  to  be  best. 

.  .  ’  *  * 

-*  6.  Motor  Parameters 

Performance  of  identical  injection  systems  may  still  vary  from 
one  motor  to  another,  making  it  necessary  to  determine  the  effect  of  certain  motor 
parameters  on  TVC  performance.  These  parameters  can  be  grouped  as:  chamber 
conditions,  propellant  properties,  and  nozzle  size  and  shape. 

The  performance  of  inert  liquids  has  been  found  to  be  quire 

sensitive  to  chamber  pressure,  as  was  indicated  in  Figure  15C ,  but  with  some 

30 

reactive  injectants,  little  effect  was  noticed.-'  A,  good  example  of  the  larger 

variation  in  performance  that  can  be  experienced  with  inert  injectants  as  the 

/  ' 

chamber  pressure  changes  is  shown  in  Figure  24,  where’  Freon  U4b„  ¥as  injected 

35  ^  t 

at  a  constant  flow  rate  during  motor  tail-off.  As  the  chamber  pressure  dropped 
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from  4-5Q  to  250  psia,  the  side-fore;  ratio  increased  by  a  factor  of  3  and  the  * 
injectant  specific  impulse  by  a  factor  of  2 ,  At  this  point,  flow  separation 
.occurred  at -the  nozzle  exit,  and  then  moved  upstream  into  the,  nozzle  as 'the  .  , 
chamber -pressure  further  .decreased .  This  separation  caused  a  sharp  drop  in  peri" 
formance  as  the  effective  nozzle  wall  area  covered  by  the  pressure  disturbance 
was  reduced^  -  , 

-h  ,Tiie  increase  in  performance  with  a  decrease  tin  chamber  pressure 
can  most  likely  be  attributed  to  greater  penetration  of  the  liquid.-  jet ,  increased 
vaporization,  and  greater  expansion  of  vapor,  resulting  in  greater  initial  flow  u 
deflect ion, and  higher  downstream  pressure  differences.  Since  the  inert  liquids- 
appear  to  be  more  dependent  on  chamber  pressure  than  certain  reactive  injeetants, . 
the  vaporization  effects  must  be  the  dominant  ones’.  '  .-y  ..  r  '  . . 

»  *  4  ‘  -  '  ' 

,  '  *  '  f 

-  -  Higher  chamber  temperati  res  would  be  expected  to  increase  the, 

performance,  of  all  liquid  injeetants,  since  both  the  vaporization  and -'chemical 
reaction  processes. are  temperature  dependent  to  a  large  extent.  Reactive' 
injectants  .will,.,  in  addition,,  be  strongly  dependent  on  the  propellant  composition. 
For  instance,,.. a  propellant  producing  large  amounts'  of  and  CO  will 'contribute 
greatly  to  the  effectiveness  of  an  oxidizer-type  injectant.  -  ■  •  ’  . 

»  «  -  * 

Size  and  shape  of  the  nozzle  is  an  important  factor  to  be  consid¬ 
ered  since  the  geometrical  effects  play  such  a  large  role  in  the  side  force 
production.  The  dependence,  of  the  optimum  injector  location  bn  the 'nozzle  L/D 
has  already  been  discussed..  /.  >-  '/* 

"  .  .  I  ' 

*  '  s  .  - 

Another  significant .parameter  is  the  nozzle  expansion  ratio . 

figure  25  compares  the  performance_at  various  injector  locations  of  a  30:1  * 

contoured  nozzle  with  another  nozzle  having  the  same  contour,  but  cut  off  at  a  -  - 

37  *’•*'  ’ 

12:1  expansion  ratio.  At  a  given  location  (  e  .).,  ,the  nozzle  with  the  greater"' 

»  1  •  -  1  ’  '  .' 
expansion  ratio  will  generally  produce  the  greater  side  force  because  of  the  'S 

greater  length  from  injector  to  exit,  and  thus  greater  available  wall  area  for  the" 
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pressure  disturbance .  However,  the  larger  nozzle  will  also  be  more  subject  to  Shock 
reflection  and  curvature  losses,  which  would  make  the  shorter  nozzle  more  effective 

a,t ..locations  close . to  the  throat.  At  the  same  relative  location  (x/L),  the  per- 

< "  *  .  "  •-  -  •* 

formance  of  . the  smaller' expansion  ratio  nozzle  will  be* greater  because  of  smaller 
Mach  numbers,  and  thus  larger  pressure  differences  and  steeper  shocks.  Both,  of 
these  trends  are  substantiated  by  the  data  plotted  in  Figure  25 •  > 


s.  The  actual  nozzle  contour  also  appears  to  have  a  substantial  effect 

on  TVO.  performance;,  Figure  26c  shows  that  at  the  optimum  injector  location,  consid- 
erably  higher  injectant  specific  impulses  for  both  inert:  and  reactive  liquids  were 

,  ,  a 

produced  in  a  contoured  nozzle  as  compared  with  conical  nozzles .  A  possible  reason 
for  this  may  be  the  greater  primary  flow  mass  concentration  near  the  wall  in  a 
contoured .nozzle,  which  would  allow  a  greater  effective  flow , deflection  by  an 

obstruction- on  the  wall.  It  may  be  possible  to  actually  optimize  a  nozzle  contour 

■  —  29 

for  TVC  performance  and  with  such  an  objective,  NESCO .  has  worked  a  great  deal 

with  an  adverse  pressure  gradient  nozzle.  However,  there  is  usually  a  loss  in 

*  >  .  *  ’  , 

axial  performance  associated  with  such  a  change  in  contour,  which  in  most  'cases 
cannot  be  tolerated.  ' 

i  ’  ’ 

C.-  AXIAL  THRUST  AUGMENTATION 


’•H- 
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A  favorable  by-product  of  the  secondary  injection  concept  is  an  increase 
in  axial  thrust  during  injection.  Because  of  the  diverging  shape  of  the  nozzle  exit 

'■  -  «i 

cone,  the  pressure  force  acting  normal  to  the  nozzle  wall  also  has  an  axial  component 
acting  an  the  same  direction  as  the  thrust.  This  component  is  primarily  a  function 
of  the  nozzle  divergence  angle,  a  ,  and  does  not  depend  upon  the  azimuthal  coordinate, 
0,  as  does  the  side  force.  Thus,  the  axial  thrust  augmentation  will  continue  to 
increase  as  the  pressure  disturbance  spreads  in  the  circumferential  direction,  and 
will  be  substantial  when  the  shock  reflects  off  the -opposite  nozzle  wall.  Also, 

the  closer  the  injector  is  to'  the  throat,  the  greater  the  projected  wall  area  in 

*  (  - 

the  cross-sectional  plane,  so  the  greater  the  thrust  augmentation. 
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1*1,  C,  Axial  Thrust  Augmentation  (cont  ) 

<  ' 

The  momentum  of  the  injected  liquid  may  also  have  an  axial  component, 
depending  on  the  injection  angle,  however,  since  the  momentum  is  usually  less 
than  10%  of  the  side  force,  this  component  would  not  be  significant  in  most  cases. 

Figures  27  and  28  show  some  thrust  augmentation  data  for  Minuteman 
subscale  nozzles^^ an<i  ^  and  illustrate  the  effects  of  injectant  flow 
rate,  injector  location,  and  angle  of  injection.  As  would  be  expected  ,  the  ' 
greater  pressure  differences  and  larger  disturbance  area  associated  with  greater 
flow  rates  cause  a  corresponding  increase  in  axial  thrust  augmentation.  Unlike 
the  side  force,  however,  at  the  locations  close  to  the  throat,  the  thrust  increment 
continues  to  increase  at  the  higher  flow  rates,  instead  of  leveling  off  and  actually 
decreasing  because  of  curvature  effects.  Also,  there  is  no  optimum  injector 
location  as  with  the  side  force,  so  that  a  location  closer  to  the  throat  always 
results  in  greater  thrust  augmentation.  It  can  be  seen  in  Figure  28  that  upstream 
injection  angles  produced  greater  thrust  ir.  crements .  This  indicates  that  the 
greater  pressure  forces  more  than  compensated  for  the  negative  momentum  axial 
components  associated  with  the  upstream  angles .  However,  it  is  possible  that  at 
locations  near  the  exit  this  relation  may  be  reversed  and  a  negative  thrust 
increment  produced. 

This  thrust  augmentation  characteristic  of  secondary  injection  is  an 
important  point  to  be  considered  when  comparing  this  system  with  mechanical  methods 
of  TVC  for  which  there  is  usually  a  loss  in  axial  thrust . 

D.  SCALE  EFFECTS  AM)  DYNAMIC  RESPONSE 

A  question  that  has  yet  to  be  fully  answered  is  whether  there  is  a 
difference  in  the  IVC  performance  between  two  different  size  motors  with  geometri¬ 
cally  similar  nozzles  and  TVC  system  parameters.  At  first  glance,  it  would  seem 
that  such  a  difference  would  exist  since  such  phenomena  as  penetration,  atomization, 
t-oundstry  layer  separation,  and  chemical  reaction  would  not  appear  to  be  directly 
scalable.  However,  experimental  data  from  Aerojet,  with  Minuteman  motors,  and 
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1U}  D,  Scale  Effects  and  Dynamic  Response  (cont.) 


at  Lockheed^-’  has  failed  to  indicate  any  definite  scale  effects »  Such  effects 
may  still  exist  though  for  different  parameter  ranges,. and,  on  the  basis  of  the  above 
phenomena,  a  slight  loss  in  performance  with  the  larger  nozzles  might  be  expected. 
Additional  experimental  verification  is  needed  to  settle  this  particular  question. 

•  ..  As  c  mpared  with  mechanical  TVC  systems,  secondayr  injection  has  an 

extremely. fast  respons or  a  short  time  lag  between  command  and  application  of 
the  control  force.  Because  of  test-stand  inadequacies,  actual  measuremenb  Oj.  .this 
response’  has  been  very  difficult .  However,  a  technique  was  developed  at  Aerojet  , 
using  exit-cone  pressure-tap  readings  instead  of  load-cell  measurements •  In  tests 
with  a  sinusoidal-injectant  flow-rate  variation,  i'  was  found  that  the  phase  lag 
between  the  side  force  and  the  flow  rate,  excluding  valve  dynamics,  was  negligible 
up  to  20  cps .  Also,  the  side  force  amplitude  ratio  varied  little  from  unity  for 

this  same  range. 


It  is  hoped  that  the  data  compiled  in  th?  a  section  on  performance 
characteristics  will  lead  to  a  better  understanding  of  the  liquid-injection 
nrocess,  and  that  the  parameter  effects  illustrated  may  be  useful  for  preliminary 
design  purposes  until  an  adequate  analytical  or  empirical  performance  prediction 
method  is  developed. 
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ANALYTICAL  -APPROACHES 


'  V 


SUMMARY  OP  PREVIOUS  WORK 


Although  the  liquid  injection  concept  has  been  investigated  for  a 

f  1  f  - 

neater  cf  yea rc,  an  analytical  model  for  adequately  predicting  the  'performance 
cf  any  ILTVC  system  has  not  been  found.  This  is  not  surprising  when  one  considers 
t.*e  complicated  processes  of  Liquid- jet  penetration,  atomization,  vaporization, 
mixing,  and.,  reaction  that  must  be  analyzed,  in  addition  to  the  extremely  complex 
three-dimensional  flow  field,  however,  from  . the  models  that,  have  been  proposed’, 
a  greater  insight  into" the  various  problem  areas  has  been  gained,  which,  along , 
vd/h  the  accumulation  of  additional  experimental  data,  leads  to  the  formulation  of 
more  realistic  assumptions  in  future  approaches.  *  •  . 


Early  in  i960,  the  first  liquid-injection  analytical  study  was  begun 

1*7 

ty  NESC-0  (National  Engineering  Science  Company)  under  contract  to  Aerojet. 

The  model  proposed  in  this  study  is  illustrated  in  Figure  29.  The  liquid  was 
assumed  to  flow  along  the  nozzle  wall,  instead  cf  penetrating  into  the  main  stream, 
and  to  form  a  uniform  distribution  of  mass  sources  due  to  vaporization.  Linearized 
theory  for  supersonic  aerodynamics  was  used  to  calculate  the  side  force  produced  by. 
the  resulting  wedge-shaped  vapor  body.  This  model  produced  fair  agreement  with 
experimental  data  for  low  injectant  mass  flow  rates,-  when  penetration  was  relatively 
small,  but  the  correlation  was  quite  poor  through  the  remainder  of  the  mass  flow 
range .  ~ 


A  much  more  sophisticated  approach  to  the  liquid  vaporization  process 

nO 

was  presented  by  NESCO  in  1961  under  a  subsequent  Aerojet  contract .  In  this 

model,  shown  in  Figure  30>  the  intermediate  processes  of  atomization  and  droplet' 
vaporization  were  included,  which,  along  with  the  droplet  trajectories,  determined 
the  resulting  vapor  "body".  The  liquid  jet  was  assumed  to  break  up  immediately  on 


entering  the  nozzle  and  a  mean  droplet  size  was  determined  from  the  empirical  rels- 

12  '  >  •* 

tions  of  Ingebo  and  Foster,  based  on  free -stream  conditions  behind  a  normal  shock 
jast  upstream  of  the  injector.  The.  droplet  vaporization  rate  was  taken  from  Penner,^ 
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-IV j  A,,  Summary  of  Previous  Work  (cont.) 
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and,  by  considering  the  drag  on  the  mean  size  droplets  due  to  the  primary  flaw,  the 

trajectory  of  these  droplets  was  computed.  The  free-stream  conditions  governing 

the  droplet  motion  and  vaporization  were  taken  as  those  existing  at  the  injector 

before  injection,  and  were  assumed  to  be  unaffected  by  the  presence  of  the  droplets 

and  vapor.  The  injected  liquid  was  thus  distributed  within  the  nozzle  as  a  continuous 

line  of  mean  size  droplets  along  the  calculated  trajectory.  The  vapor  from  this  line 

source  was  assumed  to  form  a  blunt-nosed  half-body  of  revolution  lying  along  the 

nozzle  wall  between  the- injector  and  exit  plane,  and  acting  as  a  solid  physical 

obstacle.  '  r 


i 


I  A 


♦ 


The  separated  region  upstream  of  the  injector  was  assumed  to  have  the 
shape  of  a  half- cone  tangent  to  the  nose  of  the  vapor  body.  The  equivalent  obstacle 
was  then  a, sharp-nosed  body  of  revolution,  with  the  axis  of  symmetry  being  the  nozzle 
element  passing  through  the  center  of  the  injector.  The  resulting  shock  and  flow 
field  were  assumed  to  be  symmetric  about  this  same  axis,  and  were  calculated  by  a 
simple  shock  expansion  method  for  uniform  and  constant  free-stream  conditions.  The 
side  force  was  then  calculated  by  integrating  the  wall  pressures  behind  the  shock 
and  adding  the  momentum  component  of  the  liquid  jet .  .  The  effects  of  nozzle  curvature 
in  the  cross-sectional  plane  were  not  included,  although  the  nozzle  divergence  angle 
was  taken  into  consideration. 


The  side  forces  calculated  in  this  manner  were  generally  30  to  5 0$ 
below  experimental  values,  but  this  analysis  contributed  a  great  deal  to  the  state 
of  the  art  by  indicating  the  basic  processes  that  need  be  considered,  and  by  repre¬ 
senting  the  actual  flow  field  within^the  nozzle  with  an  equivalent  simplified  model. 

20 

This  same  model  was  later  improved  at  JPL  by  including  a  more  realistic 
vaporization,  rate  and  a  relatively  sophistical0-’  droplet  drag  coefficient.  Also, 
the  separation  zone  was  modified  to  more  closely  represent  the  observed  double-shock 
pattern  as  shorn  in  Figure  2a.  Calculated  results  were  compared  with  experimental 
values  for  only  two  cases,  being  4$  high  in  one  case,  and  20$  high  in  the  other. 


TM-16-SR0 
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Broadwell^  in  1962  proposed  an  entirely  different  mod^L  for  secondary- 
,  injection  TVC,  applicable  to  both  liquid  and  gas  injectants.  This  model  is  based 
on  the  Dlast-wave  analogy,  which  transforms  results  derived  from  one-dimensional, 

,,  *  unsteady  flow  to  two-dimensional,  steady  flow.  Thus  empirical  relations  found  for 
'  an  expanding  cylindrical  shock  wave  produced  by  the  explosion  of  a  line  charge  are 
applied  to  the  problem  of  determining  the  steady- state  axisymmetric  shock  about  a 
blunt  body  of  revolution.  For  a  solid  body,  the  drag  forces  on  the  body  are  equated 
to  the  energy  released  by  the  original  line  charge  to  determine  the  size  and  shape 
of  the  shock.  For  secondary  injection,  the  injectant  momentum,  along  with  energy 

0 

eitner  taken  from  the  free  stream  to  vaporize  the  liquid,  or  added  to  the  free 
stream  due  to  a  chemical  reaction,  must  also  be  included. ■ 

Broadwell  considered  only  the  first-order  blast-wave  solution,  which 

m 

is  limited  to  very  high  free-stream  Mach  numbers.  A  second-order  solution  has  been 
.  used  by  Lahm^0  at  Vidya  in  an  effort  to  improve  the  treatment.- 

* 

Unfortunately,  the  blast-wave  analogy  is  strictly  applicable  only  for 
constant  free-stream  conditions  in  front  of  the  shock.  Also,  since  vaporization  is 
assumed  to  ta Ke  place  immediately  as  the  injectant  enters'  the  nozzle,  the  effects  of 
“  the  parameters  which  influence  the  atomization  and  vaporization  processes  (i.e.,  port 
diameter,  angle  of  injection,  and  injectant  properties),  as  well  as  the  effects  of 
boundary  layer  separation  and  nozzle  wall  curvature,  are  not  included.  Consequently, 
the  calculated  results  agree  well  with  data,  at  the  low  injectant  flow  rates,  but  the 
analysis  does  not  predict  the  decrease  in  side  specific  impulse  with  increasing  flow 

% 

rate,  nor  the  correct  variation  of  side  force  with  injector  location. 

At  Aerojet,  the  interest  in  a  theoretical  liquid- injection  model  con- 

18 

tinned  after  the  initial  investigations  by  NESCO  were  completed.  Using  the  same 
generaz  approach  as  UESCO,  modifications,  were  first  made  to  the  vapor-body  calcula¬ 
tions  no  improve  the  portions  dealing  with  the  vaporization  rates  and  drag  coeffi- 

20 

.  laese  modifications  were  very  similar  to  those  later  published  by  JPL. 
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IV,  A,  Summary  of  Previous  Work  (cont.) 

Also,  a  method  was  sought  for  simplifying  the  tedious,  time-consuming  flow  field 
calculations  included  in  the  NESCO  report.  A  dimensional  analysis  approach  was 
first  taken  in  which  the  side  force  was  related  to  the  size  of  the  vapor  body  and 
several  nozzle  and  mouor  parameters.  With  the  vapor  body  assumed  to  form  a  hali- 
body  of  revolution  for  a  single-port  injector,  the  vapor  body  radius  at  the  nozzle 
exit  was  taken  as  the  significant  body  dimension.  This  radius  was  a  function  of 
such  variables  as  the  injectant  properties,  port  diameter,  angle  of  injection, 
free-stream  properties,  and  nozzle-wall  length  between  the  injector  and  the  exit 
plane,  and  was  calculated  with  a  computer  program  similar  to  the  one  described  in 

Appendix  A.  Using  dimensional  analysis,  the  following  expression  was  derived  for 

'  22 
a  single-port  injector  and  nonreactive  injectant. 


where 


l  F 


COS  CU  /  Vs 

cF 


=  the  induced  side  force,  or  the  side  force  caused  by 

interference  only  (not  including  the  injectant  momentum); 

=  axial  thrust; 

=  nozzle  expansion  ratio; 

=  expansion  ratio  at  the  injector; 

■w 

=  nozzle  divergence  angle; 

=  nozzle  thrust  coefficient; 

=  vapor  body  radius  at  nozzle  exit; 

=  nozzle  radius  at  exit. 


The  total  side  force  was  found  by  adding  the  momentum  of  the  injectant 
entering  the  nozzle,  Fffi,  to  the  interference  force,  F  * .  The  momentum  can  be 
written  as 


_  te 

t 

\ 


+pi  ^ 


cos  9 
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IV,  A,  Summary, of.  Previous  Work  (cont.) 


where 


,=  the  injectant  mass  flow  rate;  ,  •"  - v  .  ■"  '-*■ 

=  the  velocity  of.  the  injectant  entering ‘the  nozzle;  - 

=  the  jingle  of'  injection  with  respect’  to ’the ,  normal  to'  the 
nozzle  center  line;  ,  -•  - 

V  '  .  4  *■ 

=  the  injectant  pressure  at  the  poit’exit;''  ' 

.  =  the  port#  area.  t 


The  side  forces,  calculated  in  this  manner  ‘correlated’ fairly  well  'with*  . , '  ; 

>•  ’  '  -1  ’  f  ■  '  „  f  . V  »  „  '  -  .  'V  ’ 

a  wide  range  of  experimental  data  as  shown  in  Figure  31.  "This  correlation, ‘was  ;  (  ■ 

;  ■  ■  ••  -f. 

’  encouraging  since  it  indicated  the  ability  of  the  vapor -body  model  to  predict  the  .  '•  ‘ 

*  -  '  -  i  -•  r  .  .  ■;*  '  /;  •  , 

relative  changes. produced  by  a  variation  in  injection  or.  nozzle  parameters.  '-’  • 

'M  ■  V  :  .  ■ 

'  .  •  ••  ' 

,  - .  -  1 *-  ’ 

This  same  approach  was  used  in  Reference  2-3,  to  include  the ’effects,  of 

reactive  injectants  and  multiport  injectors.  The  vapor,  body' radius,  however,  was. 

*  *  '  > 
calculated  in  a  different  manner.  An.  initial  mixture  of  yapor  and  free  stream  gases 

,  ’  ,  V, 

*  ,  ,  -  t 

was  determined  by  th ;  penetration  of  the  liquid  jet;,  and  the^  volume,’  increase  caused 

*  *  ,  *  '  ”  r  *  ,  .  >  <*  f  \  "  e 

by  chemical  reaction  was  calculated  by.  assuming  the  mixture  ,to  attain  thermochemical 
equilibrium.-  This  model,  therefore,  also  includes  the  effect  of  the  propellant 
properties,  as  well  as  additional  injectant  properties,"  and  the  circumferential 
extent  of  the  injector.  Figure  32C  shows  that  the  correlation  with  experimental  data 
obtained  with- this. model  is  similar  to  that  found  with  the  ’previous  model,,  which  was 
applicable. to  a  single  port  and^  nonreactive  injectant.  The  results  from  both 'models, 
however,  were' not  sufficiently  accurate  for  design  purposes,  especially  in  the  high 

t  '  _ 

mass  flow  range  and  for  injector  locations  far  from  the  exit.  '  As  a  consequence,-.'.  _  - 

continued  efforts  at  Aerojet  have, been  directed  'towards  the.  development  of  an  improved 

^  \  1  t 

analytical  model  for  better  parametric-  correlations .  —  .  >  ' 
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B,  CURRENT  AEROJET  PRELIMINARY  MODEL 


1.  ‘  Description 


Because  of  the  complexity  of  the  interference  pattern  created  by 
the  injection  of  a  secondary  fluid  into  a  nozzle,  a  detailed  flow-field  solution 
appears  to  be  the  only  way  to  properly  evaluate  the  effects  of  the  large  number 
of  significant  parameters  over  the  effective  range  of  each.  .-Such  a  solution  would 
have  to  be  more  accurate  and  less  cumbersome  than  the  ones  proposed  in  References  18 
and  20. 

With  this  in  mind,  a  method  of  characteristics  computer  program 
was  developed  for  calculating  a  two-dimensional  (plane)  or  axisymmetric  curved  shock 
wave  and  the  rotational  flow  behind  the  shock.  This  program,  which  is  discussed  in 
more  detail  in  Appendix  B,  requires,  as  input,  initial  conditions  along  a  start  line 
in  the  undisturbed  region  upstream  of  the  shock,  and  the  shape  of  the  solid  body 
producing  the  shock.  If  the  shock  is  detached,  an  additional  start  line  must  be 
specified  in  the  supersonic  region  between  the  shock  and  body,  with  the  subsonic  flow 
upstream' of  that  point  calculated  independently.  The  start  line  in  the  undisturbed 
region  enables  the  program  to  handle  nonuniform  free-stream  conditions  before  the 
shock,  and  thus  the  actual  initial  conditions  existing  in  a  nozzle  can  be  taken  into 
account . 


If  the  shock  strength  does  vary  along  the  shock,  the  entropy  rise, 
and  thus  total  pressure  drop,  across  the  shock  will  also  yary.  Therefore,  the 
equations  for  rotational  flow,  in  which  the  entropy  is  constant  along  a  streamline, 
but  may  vary  from  streamline  to  streamline,  must  be  used  in  calculating  the  flow 
field  betweexi  the  shock  and  body.  The  general  logic  used  in  propagating  the  solution 
is  described  in  Appendix  B. 
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li'„  Cerent  Aerojet  Preliminary  Model  (con4; . ) 

The  present  preliminary  analytical  model  for  liquid-injection 
I7C  consists  primarily  of  three  parts: 

a.  vapor  body  calculation; 

b.  flow  field  solution; 

c.  pressure  integration  over  nozzle  walls. 

< 

The  vapor-body  calculation  is  based  on  the  approach  previously  described  in 

18 

connection  with  NESCO's  model*  but  includes  the  improvements  in  vaporization 
rate  and  drag  coefficients  described  in  Reference  22.  An  additional  modification 
nas  been  made  to  include  variation  in  free-stream  properties  in  the. axial  direction, 
lo  oriefly  review  this  approach,  a  mean-size  droplet  is  found  from  empirical  rela¬ 
tions,  and  the  'trajectory  of  this  droplet  is  calculated  by  considering  the  drag 
forces  acting  on  the  droplet  because  of  the  free  stream.  The  presence  of  the 
droplets  is  assumed  to  have  no  effect  on  the  free-stream  properties  used  to  cal¬ 
culate  the  drag  forces.  As  the  droplets  travel  along  their  trajectory,  vaporiza¬ 
tion  occurs  because  of  the  high  free-stream  temperatures.  The  resulting  vapor 
body  is  found  by  considering  the  injectant  distributed  along  the  mean  size  droplet 
trajectory  to  act  as  a  line  source,  and  by  assuming  the  vapor  to  form  a  half  body 
of  revolution.  This  entire  procedure,  described  in  more  detail  in  Appendix  A 
along  with  the  derivation  of  the  significant  equations,  has  been  incorporated  into 
a  computer  program  on  the  GE  225,  In  Figure  33,  calculated  vapor  bodies  are 
compared  with  erosion  patterns  formed  by  injection  of  Freon  ll4B2  into  a  Minuteman 
no?  >le  during  a  hot  firing.  If  the  actual  vapor  body  is  axisymmetric,  these  com¬ 
parisons  indicate  its  size  and  shape  are  represented  fairly  well  by  the  analytical 
model. 


Before  the  vapor-body  shape  can  be  input  to  the  flow  field  solu¬ 
tion,  it  must  be  modified  to  include  the  separated  region  upstream  of  the  port. 

For  the  saKe  of  simplicity  in  this  preliminary  model,  this  region  was  first  assumed 
to  be  conical  and  tangent  to  the  nose  of  the  vapor  body  as  in  Reference  18,  Thus 
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only  e.  single,  attached  shock  results,  as  shown  in  Figure  3a,  and  the  entire  flow- 
field  can  be  calculated  with  the  method  of  characteristics  program  previously 
described.  Provision  can  be  made  later  to  include  the  general  two  shock  system 
shown  in  Figure  2a,  or  the  limiting  bow  shock  at  the  high  mass  flows  as  shown  in 
Figure  3b . 

To  find  the  vertex  angle  of  this  conical  separation  region,  the 
24 

relation  derived  by  Mager  was  used  in  which  the  plateau  pressure  rise  associated 
with  the  separation  of  a  turbulent  boundary  layer,  I^/lco,  given  as  a  function  of 
the  free-stream  Mach  number,  I^so,  the  flow  deflection  angle,  6,  and  the  ratio  of 
specific  heats,  y.  The  resulting  equation  is 


where  K  is  the  ratio  of  the  square  of  the  Mach  number  at  the  separation  point  to 
2 

Moo  *  Experimental  data  indicates  this  ratio  to  be  approximately  0.55 . 


For  given  values  of  Moo  &«d  y,  the  pressure  rise  across  a  conical 
shock  can  be  found  as  a  function  of  cone  angle.  By  satisfying  this  additional 
relation  along  with  the  above  equation,  the  resulting  value  of  p^/p^and  6  can  be 
found.  Figure  34  shows  the  variation  of  these  variables  with  M^for  y  =  1.2.  The 
value  of  6  determines  the  size  of  the  separation  region  since  it  is  assumed  to  he 
tangent  to  the  vapor  body. 


To  utilize  the  method  of  characteristics  flow  field  solution, 
the  assumption  is  made  that  the  shock  and  flow  behind  it  is  axisymmetric  with 
respect  to  the  nozzle  element  passing  through  the  center  of  the  injector.  This 
is  probably  the  weakest  point  in  the  present  model,  since  the  actual  flow  is 
three-dimensional.  However,  it  was  felt  that  the  axisymmetric  solution  should 
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represent  the  actual  flow  to  a  fairly  good  degree,  and  could  "be  corrected  for  the 

t  major  three-dimensional  effects.  A  typical  calculated  vapor  body,  separated .region, 
shock,  and  flow- field  network  is  shown  in  Figure  35.  '  '  ' 

0  «, 

The  last  step  in  finding  the  side  force  is  the  integration  of  the 
calculated  pressures  over  the  nozzle  trails .  By  superimposing  the  given  nozzle 
shape  upon  the  axisymmetric  flow-field  solution,  the  shock  and  nodal  points  in  the 
characteristic  solution  can  be  located  along  the  nozzle  wall.  Figure  36  shows  a 
.ypical  intersection  of  a  calculated  axisymmetric  shock  wave  with  a  nozzle  wall, 
lo  calculate  the  side  force,  these  points  are  then  projected  into  a  plane  parallel 
to  the  nozzie  center  line  and  perpendicular  to  the  plane  formed  by  the  nozzle  center 
line  and  the  shock  axis  of  symmetry.  Thus,  when  the  calculated  pressures  are  inte¬ 
grated  in  this  plane,  the  effects  of  the  nozzle  divergence  angle  and  wall  curvature 
in  the  cross-sectional  plane  are  included.  Also,  provision  has  been  made  for 
handling  the  case  of  a  shock  intersection  with  the  nozzle  wall  above  the  center  line, 
which  produces  a  negative  force  component. 

) 

The  axial  thrust  augmentation  is  found  in  the  same  manner,  except 
the  solution  points  are  projected  into  a  plane  normal  to  the  nozzle  center  line. 

In  this  case,  no  reversal  of  force  occurs  when  the  shock  crosses  over  the  center 
line. 

An  additional  item  calculated  is  the  center  of  force,  or  the 

effective  point  of  application  for  the  induced  side  force,  F  This  axial  distance, 

s 

xc,  (using  the  nozzle  throat  plane  as  the  reference  plane)  can  be  found  from  the 
following  expression: 

p  As 

I  x  (p-pco  ) 

X  = 

^  As 

(P-Poo )  *A 


As 


u  o 


x  (P-PTO  )  dA 


(4) 
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•As  =  ..the  nozzle  wall  area  behind  the  shock  in  the  jalane 
£  -«  perpendicular to  the  side  force; , < .  . 

x  =  the  axial  distance  of  dA  from  the  reference  plane; 

P  ='  the  average  static  pressure  acting  upon  dA  during 
■  v .  injection;  .  •  j  , 

=  'the. ‘average  staoic  pressure  that' would  act- on.  dA 
with  no  injection. 
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■  This  pressure  integration  routine  has  been.’ combined 'With  the  flow  ' 

■  '■  ■  '  s’  V.-'-  ‘  :  '  ‘  -  • 

field-,  solution  into  a  single,  computer  program,  which  takes  less  than  1  minute'  per 
rim  on  the.  IEM  7094,  The.  calculation  of  the  vapor  body  and  separation  ..region  will  , 
also  be  included  in -this  one,  program  when  the  final  form  of*  the.  body ’shape  has  been 
decided.-  < 
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To  evaluate  the  performance  prediction  capability  of  this  model, 

Ni  *  /  ,  .  ’  T  * 

side  forces,  were  calculated  for  a  number  of  different  motors  encompassing  a  wide 

'  “  .  •  „  *  J.  .  ;  *  . 

range  of  injection  and  nozzle  parameters.  Axial  thrusts  for  these  motors  varied 
from  l,00Q.4to  100,000  lb,  ,  nozzle  expansion- ratios  from  6.64:1 ( to  24.8:1,  ftnjeetant 
mass  flows  'for  Freon  113  .and  114b1  from  2$.to‘22 $  of*  the  primary  'flow,  and  injector 
locations  at  expansion  ratios  from  3:1  bo  12.5:1.  When  these  .calculated  results  . 

^  ,  *.  1  '  '  *  ,  /  %  -*  if,* 

were  compared' .with  the -measured  data  for  these  motors,  it  was  found  that  the  model 
predicted  side  forces  that  were  generally  40  to  50$  too  low.  -  However,  the ’ratio  v 
of  calculated  to  measured  side  forces  did  not  appear  affected'  by  the  wide  parameter 


r*;i 


- 

v*jJ 


\  ' 


variation. 


>  i  The -main  reason  for  these,  low' values  .can  be  found  by  comparing1* 

~r  .  .  z  .<  s  .  ‘  ‘  ‘  ’  -  •.  *  ,  -M. '  4  .  *’ 

the  calculated  and  measured'pressrure  distributions .  Figure ,37 'shows  that  the  >' 
theoretical  axial  pressure. di str ibut  ions  along  the  vapor  body  itself  agree  fairly’  *: 
well  with  the  measured  values,  but  the  circumferential  spread  of  the  calculated  1 

disturbance  is  much  less  than  indicated  by  experimental  data. '  This  latter  result  ' 

<  -  ■.  * 
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IV;,  B,  Current  Aerojet  Preliminary  Model  ( cont . ) 

is  more  evident  in  Figure  38?  where  the  traces  of  both  the  theoretical  and  experi- 
,  •  mental  shock  raves  along  the  nozzle  wall  are  shown.  The  largest  deviation  between 

the  two  shock  patterns  occurs  in  the  separation  region,  as  might  be  expected,  with 
*  tne  measured  disturbance  hyperbolic  shaped,  while  the  calculated  shock  and  separa¬ 
tion  form  a  vortex  of  a  triangle.  The  assumed  conical  shape  for  the  separation 
region  is  thus  obviously  incorrect.  Farther  downstream,  the  difference  might  be 
attributeu  to  three-dimensional  effects  and/or  too  small  a  calculated  vapor  body. 

In  addition,  boundary  layer  separation  may  also  occur  in  the  area  downstream  of  the 
port  due  to  the  impingement  of  the  shock  on  the  nozzle  wall.  This  would  tend  to 
spread  the  disturbance  farther  upstream  from  the  shock  itself.  Since  relatively 
*  high  pressures  occur  in  this  additional  area  not  covered  by  the  theoretical  shock, 
the  difference  in  area  causes  a  substantial  difference  in  side  force.  This  pre- 
„  liminary  model  has  thus  Indicated  definite  areas  that  must  be  improved  before  the 
final  form  of  the  model  is  determined.  - 


V. 


RECOMMENDATIONS  FOR  FUTURE  WORK 


Although,  the  results  from  ihe  preliminary  model  are  not  satisfactory  at 
present,  the  consistency  of  the  relation  between  the  calculated  and  measured 
side  forces  is  very  encouraging.  This  model  thus  not -only  forms  a  firm  basis  for 
future  work,  but  also  indicates  the  specific  areas  needing  improvement  and  the 
definite  directions  to  go  in  each  area.-  It  is  estimated  that  the  model  in  its  present 
form  is  approximately  70$  complete.  To  bring  it  to  full  completion,  the  following 
tasks  will  be  undertaken: 

1.  a  better  representation  of  the  separation  region; 

2.  corrections  to  the  flow  field  solution  for  three-dimensional 
.  effects; 

3.  more  realistic  vapor-bo<?y  calculations. 

These  tasks  will  depend  greatly  on  basic  experimental  data  in  the  form  of 
pressure  distributions  and  visual  evidence,  such  as  shadowgraphs  and/or  schlieren  photo¬ 
graphs  of  the  flow.  A  current  IR&D  program  for  gas-injection  TVC  has  been  designed 
to  provide  specific  information  of  this  kind,  which  will  be  directly  applicable  to 
the  first  two  tasks. 

According  to  data  now  available,  the  separation  region  apparently  spreads 

out  along  the  wall  much  more  than  would  be  1  licated  ,by  a  simple  half -conical 

shape.  Thus,  the  flow  over  this  region  is  three-dimensional,  and  cannot  be  included 

in  the  axisymmetric  flow  field  solution.  Therefore,  it  will  be  necessary  to  calculate 

this  region  separately,  possibly  by  semiempiricsl  means.  In  addition,  in  the  higher 

inject  ant  flow  ranges,  when  a  strong  shock  is  produced,  the  oalc’ilation  of  this  shock 

and  the  subsonic  flow  behind  it  must  also  be  included.  If  this  subsonic  region  is 

25 

relatively  small,  a  simplified  approach  such  as  the  method  of  Moeckel  might  be 
adequate  for  defining  the  initial  conditions  for  the  supersonic  solution.  If  a  more 

26  i 

sophisticated  subsonic  solution  is  required,  the  method  of  integral  relations  has 
proven  successful  for  calculating  detached  shocks  with  the  accompanying  subsonic  flow,  7- 


t 
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and  can  be  combined  with  the  method  of  characteristics  to  give  the  entire  flow  field. 
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The  .two  major  three-dimensional  effects  for  whickthe  axisynimetric  ilaM-field  <,  V 
,  .1  at  ion  should  he  corrected  are  the  initial  conditions  in  the  nozzle  'before*  the 
snccK,  which  are  not  symmetric  with  respect  to  the  shock  axis,,  arid  the  nozzle  wall';  •  i 
•  curvature  in  the  cross-sectional  plane.  In  this  particular  area,  the  data  obtained 
from  the  aforementioned  gas-injectica  program  will  be’  quite  valuable  since  a  given 
set  of  injection  conditions  will  be  kept  constant  for  injection:  .  .  ‘v  ■* 


\ 

C 


-tic 

■  * 


(  f 


1;^.  through  a  flat  ,  plate  into  a  uniform  free- 3 1 ream;  >.■ 
2;  into  a  two-dimensional  (plane)  nozzle; 

3.  into  an  axisymmetric'  nozzle.  ’  ,  V‘  ’  vt 


Ey  comparing  the  shock  shapes  and  pressure  distributions  obtained  for  each  case, , 


the  three-dimensional  effects  will  "be  determined. 


The  ultimate  'answer  to  the  problems  .in  this,  area,  is  the  use  of  the  three- 
dimensional  method  of  characteristics  in  the  flow-field,  solution*-*  „  The  logic  involved 
in  programing  such  a  solution  is  much  more  involved  in  comparison  with  the  two- 


iimensional  or  axisymmetric  method  of  characteristics,  but  the  approach  presented. 
-28  ' 

by  Morel ti  has  been  successfully  orogramed.  Since  the  three-dimensional  subsonic 


solution  can  also  be  included,  as  shown  in  Reference  27,  the  entire  flow  field  about 

,ne  body  with  its  attached  three-dimensional  separation  region  could  once  again.be 

calculated  in  one  step.  Also’,  the  vapor  body  would  no  longer  be  restricted  to  a. 

half -body  of  revolution,  and  could  thus  include  the  effects  of  liquid- jet  penetration 

* 

ani  multiport  injectors.  Because  of  the  great  value  of  a  three-dimensional  method  of 
characteristics  program,  not  only  for  a  secondary  injection  model,  but  ror  many  v 
oti-.r  supersonic  flow  problems  as  well,  work  has  been  initiated  here  for  the  develop- 
ment  of  such  a  program.  ‘  «.  ^  1  •  • . 


The  vapor  body  calculation,,  the. heart  of  the  liquid,  injection, model,  has  , 

*  f,  t  w  . 

or.?  tig  weakness-- the  use  of  the  free-stream  conditions  outside  the  body  to  determine 

+  -  .  *  ’  s,  ■ 


droplet  trajectories,  vaporization,  etc.  In  reality,  the  droplets  are  traveling  in 
a  medium  which  is  a  mixture  of  yapor  and  free-stream  gases,.  In  a  model  proposed* by  . 
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29  - 

Mayer,  the  droplet  behavior  is  based  on  the  properties  of  such  a  mixture,  and  the 

resulting  droplet  velocities  are  much  lower  than  those  calcula+ed  by  the  present 

method.  For  a  given  amount  of  vapor  generated,  a  larger  vapor  body  will  thus  be 

produced,  resulting  in  a  greater  side  force.  The  present  vapor  body  calculation  will 

therefore  be  modified  to  include  this  more  realistic  treatment  of  the  interaction 

between  the  liquid  droplets,  vapor,  and  free-stream  gases-. 

Once  satisfactory  results  have  been  obtained  for  injection  of  an  inert 
liquid  through  a  single  port,  the  model  will  be  expanded  to  include  reactive  liquids 
and  multi -port  injectors.  The  chemical  reaction  between  the  injectant  and  free 
stream  can  probably  best  be  handled  by  thermoebemieal  calculations  similar  to  those 
described  in  Reference  23.  With  the  three-dimensional  method  of  characteristics, 
the  vapor-body  shape  can  be  made  dependent  on  the  number  and  spread  of  the  injection 
ports,  as  well  as  the  liquid- jet  penetration. 

The  tasks  proposed  here,  with  the  exception  of  the  three-dimensional  method 
of  characteristics,  will  not  amount  to  major  efforts,  but,  for  the  most  part,  must 
await  the  forthcoming  experimental  data.  When  the  three-dimensional  characteristic 
solution  has  been  completed,  and  each  portion  of  the  model  satisfactorily  checked 
with  experimental  data,  the  entire  model  will  he  combined  into  a  single  computer 
program. 
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< 

Liquid  injection  may  noi >i  be  considered  a  proven  reliable  method  of  thrust- 
t  vector  control  as  a  result  of  its  performance  on  the  Polaris  and  Minuteman  test 
»  flights.  Although  LITVC  appears  best  suited  for  upper-stage  applications,  static 
*  test  results  for  the  Titan  IIIC  120-in. -dia  solid  booster  have  shown  the  performance 
of  a  system  using  22  i-njectant  to  be  entirely  adequate.  Also,  because  of  its 
last  response,  LITVC  has  been  designated  for  the  upcoming  Sprint  anti-missile  missile. 
Therefore,  this  type  of  thrust-vector  control  will  most  likely  be  used  extensively  in 
future  designs,  with  performance  increases  coming  as  the  technology  grows.  Although 
greater  performance  can  be  attained  with  gas  inject ants ,  it  will  be  some,  time  yet  before 
a  system  witn  a  hot-gas  valve  wiL1  reach  the  stage  of  reliability  current3.y  associated 
with  liquid-injection  systems. 

Despite  the  fact  that  LITVC  systems  have  proved  successful  in  many 
applications,  there  is  still  relatively  little  known  about  the  complex  flow 
rnter action  whicn  produces  the  side  force.  The  present  system  designs  have  been 
reached  only  after  extensive  and  costly  testing,  and  still  may  not  be  completely 
•-p*,rmum.  There  is  thus  a  great  need  for  a.  method  of  predicting  the  performance  of 
any  LITVC  system,  which  can  he  used  in  lieu  of  many  tests  to  determine  more  realistic 
preliminary  designs.  Such  a  method  would  t-e  used  for: 

1.  evaluating  the  performance  of  a  given  system; 

2.  optimizing  system  parameters  for  a  given  application; 

3.  extrapolating  subscale  test  data  to  full-scale  designs; 

4.  sizing  system  components, 

* 

The  simplified  analytical  models  and  empirical  methods  that  have  "been 
.  attempted  to  date  do  not  include  sufficient  detail  for  design  purposes.  Only  a 
model  which  actually  represents  the  flow  within  the  nozzle  during  injection  can 
ir.cjude  the  first  order  effects  of  all  parameters,  while  contributing  tc  a  greater 
knowledge  and  understanding  of  the  associated  phenomena. 
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VI,  Conclusions  (cont.) 


f  • 


.  ,  It  -was  towards  this  objective  that  the  development  cf  the  model  described 

-%  * 

in  this  report  was  initiated.  When  completed,  this  model  will  have  the  additional 

'  '  .V  *  > 

advantage  of  providing  detailed  pressure  distributions  for  nozzle  stressing,  heat- 
..transfer  studies,  and  injector  design.  Also,  individual  portions  of  the  solution  can 
be^ easily  modified  as  additional  information  becomes  available. 

*  * 

i,  '  *«  - 

In  the  very  near  future,  when  the  experimental  results  from  the  current  gas- 
injection  cold-flow’  i.est  program  are  available,  this  model  should  be  completed for  ’  - 
inert  in jectants  and  single  port. injectors.  Additional. modifications  to  handle 
reactive  liquids  and' multiport  injectors  do  not  appear  to  be  major  undertakings. 

The  final  objective  is  a  single  computer  program  that  will'  include,  any  inject  ant, 
gas  or  liquid,  for  any  kind  of  injector,  in  any  type  of. "nozzle. 
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NOMENCLATURE  - 

-  Injector  port  diameter 

*v 

Ag  -  Nozzle  wall  area  behind  shock  projected  into  plane  perpendicular  to 
side  force  .  ~~  ' 

-  Droplet  drag  coefficient 

(C  )v  -  Vapor  specific  heat 

C_  -  Axial  thrust  coefficients 
F 

I 

C  -  Injector  port  discharge  coefficient 

d^  -  Injector  port  diameter 

Dq  -  Initial  mean  droplet  diameter 

D  -  Nozzle  exit  diameter 

L*  -  Nozzle  throat  diameter 

fA  -  Jbdalthr“rt 

(Fa)o  -  Axial  thrust  without  injection 

Fp  -  Drag  force  on  droplets 

F 

M  -  Component  of  injectant  momentum  contributing  to  side  force 

F  -  Side  force 

s 

F£*  -  Side  force  excluding  injectant  momentun 

Ic  -  Motor  axial  specific  impulse 

(I  )  -  Injectant  side  specific  impulse 

s  s  * 

K  -  Ratio  of  Mach  number  at  separation  point,  squared,  to  free-stream 
Mach  number,  squared 

L  -  Nozzle  length 

L  -  Injectant  latent  heat  of  valorization 
v 
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•  NOMENCLATURE  (coni.) 

r  -  Total  velocity 

-  Total  droplet  velocity 

-  Initial  injectant  velocity 

-  Relative  velocity 

-  Motor  propellant  mass  flow  rate 

-  Inject ant  mass  flow  rate 

-  Weber  number- 

-  Axial  coordinate 

-  Distance  of  center  of  force  from  nozzle  tb  ’oat 
c 

(x,y,z)-  Cartesian  coordinates 
(r,0,x)-  Cylindrical  coordinates 

-  Nozzle  wall  angle  (Mach  angle  in  Appendix  B) 

-  Vaporization  rate 

-  Vaporization  rate  in  a  stagnant  atmosphere 

-  Ratio  of  specific  heats 

-  Flow  deflection  angle  across  shock 

-  Shock  angle 

-  Nozzle  expansion  ratio  at  injector 

-  Nozzle  expansion  ratio  at  exit 

-  Azimuthal  coordinate  (Flow  angle  in  Appendix  B) 

-  Vapor  thermal  conductivity 

jx  -  Absolute  viscosity 

v  -  Kinematic  viscosity 

X,  -  Distance  along  droplet  trajectory 
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Superscripts  '' 

(  )  -  Denotes  average  quantities  in  finite  difference  solution  .  ‘ 

Subscripts 

'  *  ^ 

A,  A].,  Denotes  points  in  shock  point  solution  . 

A2>  -  *  -  - 

Bl>  b2, 

C,  D,  E,  - 

F  '  . 

d  -  Pertains  to  droplets 

i  -  Denotes  conditions  at  injector  port 

I  -  Pertains  to  left-running  characteristic  •  '  -  .  ,  • 

>  r 

(  *  1  -  -  ’ 

J  -  Pertains  to  right-running  characteristic ;  ^  ^ - 

L  -  Pertains  to  the  inject  ant  in  liquid  form  " 

n  —  Denotes  point  to  he  computed  in  finite  difference  solution 
*  ,  .  ■  '  v 
n-1  -  Denotes  point  previously  computed  in  finite  difference  solution 

N  -  Denotes  conditions  behind  normal  shock 

v  '  -  .Pertains  to  .inject ant  vapor  ’■  " 

oo  -  ; Denotes  free-s bream  conditions 

o  -  Denotes  initial  conditions 

i 

1,2.3)  -  Denote  points  in  interna?  point  and  boundary  point  solutions 
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MOTORS  PARAMETERS 


Aerojet  HOTS 
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Figure  '2 
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MOTOR  PARAMETERS 

-23. 5: X  Contoured  Nozzle 
D*  =  .4.983" 

Isp  -  279 
CF  “  1,75 
Pc  ~  500  psia 
Propellant  -  Solid 


INJECTION  PARAMETERS  * 
1200  psia  - 
Single  Port 

Injectant  -  Freon  ll4B 
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LITVC  Performance  During  Motor  Tail-off 
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(Aerojet  -  Minuteman,  R. 

KOTOR  PARAMETERS 

23.5:1  Contoured  Nozzle 
D*  =  4.938" 

Isp  “  279 
CF  -  1.75 
Pc  -•  450  psia 
Solid  Pro reliant 


INJECTION  PARAMETERS 

^  =  30? 
pinj  =  15°°  Psia 
Single  Port 

Injectant  -  Freon  114b0 
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X/L  =  C.45 


Axial  Thrust  Augmentation  due  to  Liquid  Injection 


Figure  27 ' 
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Liquid  Injection  Linearized  Model 


Figure  29 
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Correlation  of  Experimental  Data  with  Approximate  Model 
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'  The  portion  of  the  present  liquid -inject ion  model  that  calculates  the  - 
equivalent  solid’ body  formed  by  the  inject ant  vapor  follows  the  general  approach 
described  by  NE8C0  in  Reference  18.  A  number  of  modifications  have  been  made, 
however,  to  include  the  effect  of  variable  free -stream  conditions,  and  to  improve 
the  vaporization  rate  and  drag  coefficient.  Further  work  is  required  to  properly 
evaluate  the  coupling  effects  between  the  free  stream,  vapor,  and  droplets,  and 
to  include  the  effects  of  liquid -jet  penetration. 

Since  the  injectant  is  assumed  to  break  up  into  droplets  immediately  on 
entering  the  nozzle,  the  first  step  in  the  present  model  is  the  calculation  of 
the  mean  size  droplet  trajectory.  By  considering  the  droplets  along  the  trajectory 
as  a  line  source  of  mass,  the  distribution  of  vapor  is  then  found  by  assuming  a 
semicircular  cross-section.  . 


To  determine  the  diameter  of  the  mean  size  droplet,  an  emnirical  relation 

12 

derived  by  Ingebo  and  Foster  for  injection  of  a  liquid  perpendicular  to  a 
subsonic  free  stream  is  used.  The  free  stream  conditions  are  taken  as  those 
behind  a  normal  shock  located  at  the  injector. 
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This  relation  was  considered  more  realistic  than  the  well-known  lluklyama -Tanas awa 
11 

expression  ,  which  was  derived  for  injection  parallel  to  the  freestreara  and  includes 
no  effect  of  the  port  diameter. 
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An  expression  for  this  vaporization  rate  is  given  in  Reference  46  for 
quasi -steady -state  evaporation  of  a  spherical  droplet  in  a  stagnant  atmosphere, 


2  Xv 
apvPx, 


(A5) 


An  empirical  correction  can  be  made  to  this  equation  to  make  it  applicabl 

46 

for  fluid  motion  with  respect  to  the  droplet, 

1/2  1/3 

p  =  (1  +  0.276  Re  Pr  )  (A6) 

To  arrange  the  trajectory  equations  into  their  final  forms,  equation  (A4) 
is  first  rewritten  as 


2rddra  =  -pat 


(AT) 


Writing  equation  (A2)  in  component  form,  and  substituting  equation  (A4) 
into  it  gives: 


du 
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.  dr, 
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(AS) 


dv 
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^  dv 
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(A9) 


Equations  (A7),  (AS),  and  (A9)  form  a  set  of  simultaneous  ordinary 
differential  equations  to  be  solved  for  the  droplet  radius,  r^,  and  the  droplet 
velocity  components  u  and  v  as  a  function  of  time,  t. 
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To  describe  the  actual  droplet  trajectory,  the  iollowing  two 'equations 
are  then  needed, 


dx 

=  udt 

ay 

=  vdt 

(A10) 

(An) 


After  determining  the  droplet  trajectory,  the  next  step  is  the  calculation 
of  the  bcdy  formed  by  the  vapor.  The  mass  of  vapor  generated  per  drop  in  an 
interval  d^  along  the  droplet  trajectory  is 


dm  =  -  47rr  dr. 


Substituting  equation  (A4)  into  the  above  expression  gives 


(A12) 


and 


dm 

dt 

dm 

d? 


2  V  PT  r.,  P 
L  d 


dm , 
dt 


dt 

d. 


2irpLrdp 


(A13) 


The  total  number  of  droplets  per  unit  time  entering  the  nozzle  is 
ws 


S  = 


h.  3 

4-  TTr  J  P_ 
3  o  L 


(Al4) 


i 


^  !  ' 


1  X  ' 

ti 

•  U : 

r 

u 

if.  i 


i  * 


Multiplying  equation  (A13)  by  equation  (Al4)  gives  the  total  mass  of  vapor 
per  unit  time  generated  in  the  interval  dS, 


dw  _  tj  _3 

d£  d$  2" 


“s  rdP 

r  3  va 
o 


(A15) 
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S 


i; 


I:  ' 


r. 

-&J 


it 


Since  the  droplet  trajectories  are  ary  close  to  the  wall,  as  shown  in 
Figure  A1  the  droplets  are  assumed  to  form  a  lin-:  source  lying  along  the  nozzle 
wall.  The  vapor  thus  forms  a  half -body  o"  revolution  also  lying  along  the  wall. 


The  radial  velocity  of  the  vapor  at  the  surface  of  the  body  due  to  the 
line  source  is 


y  _  _ dw 

v  P  7TR  dx 


and  the  slope  of  the  vapor  body  is 


^v  _  V  co 
dx  ~  Uoo 


p  TT  R 
r  V  V 


dw 

as 


(A16) 


(A17) 


RE-SCO  in  Reference  17  shows  by  a  Taylor  series  expansion  in  pv 
that  at  the  body  surface 


P  v  V  ~  P  oo  Vco 
V  v 


(AlS) 
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Therefore,  substituting  equations  (Al6)  and  (Al8)  into  (A17)  gives 


dR 


dx 


pcotfco  TTR 


dw 

TT 


(A1 9) 


and  with  equation  (A15) 


2R  XR  -  >4” - -  —  ^ 

v  v  ^  IfeoOSo  r.3  V 


dx 


(A20) 


0  'd 


Equations  (A7)  through  (All)  and  equation  (A20)  constitute  6  equations  for 
the  6  unknowns ,  u,  v,  x,  y,  r^,  and  R^ as  functions  of  the  independent  variable  t. 
Since  this  is  a  "marching  type"  problem,  the  equations  can  be  put  into  finite- 
difference  form,  which  allows  the  numerical  solution  to  advance  from  given  initial 
values  at  successive  discrete  intervals  of  time. 

As  an  example  of  the  procedure  used  for  deriving  the  finite  difference 
equations,  equation  (A8)  is  integrated  across  a  given  interval  in  the  following 
manner . 


un 


du 


(Uoo  _u)‘ 


3 

T 


Poo  ^Cp^x 

pi  ^ 


u 


n-1 


-(r,) 
v  a'n 


(rJ  . 

a  n-l 


dr. 


(A2l) 


where  the  quantities  Poo  ,  Uoo  ,  (C^)x>  and  ""fT  ar^  taken  as  constant  average 
values  in  each  interval.  The  above  equation  can  then  be  integrated  to  give 


or 


1  1 

[uoo 

"  u 

n  / 

\ 

u 

n 

=  Uoo 

- 

poo  _  u 


n-l  / 


3  PM'Vx 
TT 


Pi 


r  ch  -  x  -  ] 


(A22) 


1  1 

.  3 

900  (Vx  f(r  )  -  M  ] 

}U0°  -  un_1j 

+  T 

pl  £  n"1  J 
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The  other  finite  difference  equations  can  be  derived  in  a  similar  manner, 
and  are  listed  below 

1/2 


=  [(ra>n-l  •  P  (t- 

,  Ks  (Ol) 


n 


=  1 


n-1 


X  =  X  .  + 

n  n-1 


X 


U  +  U  - 

n  n-1 


V  +  V  - 

n  n-1 


PL  P 

(t  -  t  , ) 
v  n  n-1' 


(A23) 

Z  Oa>  ,  -  ('a>  ]  } 

L  n-1  .  n  J 


n  ^n-l 


n 


W*..  + 

L  n-x 


3 

T 


(t  -  t  ..) 
v  n  n-1 


*  P 


(A25) 


(A26) 


(O  +  W 


1/2 


UvJ  +  (vJ  -1  11  n"1  J 


poo  uoo  r  3  L(y  )  +  (Va) 

o  n 


n-1 


(A2T) 


where 


\  =  (u2  +  v2) 


1/2 


The  initial  values  used  (@t  =  C)  are 


u 


V_  sin  9 
L 


v  =  VT  cos  9 
o  L 


r  = 


yQ=  og 


do 

2 


where  do  is  calculated  by  equation  (Al) 
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.  v  fcr 

For  the  solution,  an  iteration  technique  is  used  in  which  the '’values  for' 
ffco  ,  uoo  3  }  '(c  )  ,  and  (C^)  for  the  first  pass  in  a  given  interval,  take  , 

the  form  x  y  * 


a  «  a„ 


(A29) 


After  the  first  pass  through  equations  (A22)  (A27),  a  new  value,  is  given  to 

these  terms. 


a  ,  +  a 
n-1  n 


(A30) 


where  a .  is  calculated  from  the  results  of  the  first  pass,  and  the  solu  lon  is 
repeated.  This  procedure  continues  until  successive  values  calculated  for  u^ 
satisfy  the  following  criterion 


(u) 

xrf~  1  ^  •°001 

n  i-1 


(A31) 


The  converged  values  of  the  dependent  variables,  An,  are  then  used  as  the  initial 
values,  kf  for  the  next  interval  to  continue  the  solution. 

*  Cn~i) 


This  solution  has  been  programmed  for  the  GS225  computer  and  takes 
approximately  5  minutes  per  run.  This  program  requires  the  following  input: 

1.  Propellant  Properties 

a)  Y  ,  ratio  of  specific  heats 

b)  R,  Gas  Constant,  (ft-lbf/lbm-°R) 

c)  Pc,  Chamber  Pressure,  (psia) 

d)  Tc,  Chamber  Temperature,  (°R) 
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Injectant  Props -ties 

a)  P  T,  Liquid  Density,  (lbm/ft^) 

ii  2 

b)  Vy  Liquid  Kinematic  viscosity,  (ft  /sec) 

c)  (T  y  Liquid  Surface  Tension,  (ibm/ft) 

d)  X  ,  Vapor  Thermal  Conductivity,  (BTU/sec-ft-°R) 

e)  .  L^,  Latent  Heat  of  Vaporization,  (BTU/lbm) 

f)  T j.  Boiling  Temperature,  (°R) 

g)  Vapor  Specific  Heat,  (BTU/lbm-°R) 

Injector  Parameters 

a)  wg,  injectant  mass  flow  rate,  (ibm/sec.) 

b)  9  ,  injection  angle  (degrees) 

c)  d..  Port  diameter,  (in.) 

d)  P^j,  Injection  pressure,  (psia) 

e)  Cw,  Port  Discharge  Coefficient 

Table  of  Free  Stream  Mach  Numbers  vs.  X. 
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PLOW-FIELD  SOLUTION 
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The  method  of  characteristics  has  long  been  recognized  as  a  valuable  method 

for  solving  hyperbolic  partial  differential  equations,  and  has  been  used  extensively 

for  compressible,  inviscid,  supersonic  flow  field  solutions.  Mathematically,  this 

method  transforms  parcial  differential  equations  that  describe  the  flow  over  the 

entire  flow  field  to  ordinary  differential  equations  which  apply  only  along 

characteristic  lines  in  two-dimensional  flow,  or  characteristic  surfaces  in  three- 

dimensional  flow.  The  general  theory  of  the  method  of  characteristics  and  the 

actual  derivation  of  these  ordinary  differential  equations,  sometimes  called 

compatibility  equations,  will  not  be  discussed  here,  but  can  be  found  in  a  number 

47  48  49 

of  texts,  such  as  Courant  and  Friedrichs  ,  Ferri  ,  and  Shapiro  . 

Since  the  flow  field  for  LITVC,  as  shown  in  Figure  35,  will  include  curved 
shocks,  the  flow  behind  the  shocks  will  be  rotational,  i.e.,  the  entropy  will  be 
constant  along  a  given  streamline,  but  will  vary  from  one  streamline  to  the  other. 
From  Shapiro,  for  rotational  axisymmetric  flow,  the  characteristic  lines  are  the 
Mach  lines,  described  by  the  equations, 

m  = 

IdxJ  j 


r 


x 


tan  (  9  +  a) 

(Bl) 

tan  (0  -  a) 

(B2) 

J  Characteristic 
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and  the  compatibility  equations  applicable  along  these  lines  are 

(dr)_ 

o'?  ft  r*nc  CL  * 


cot  q  (dV\  +  sin_G 
■»r  \av/T  oTnTfl  +rr  1 


slne.Ol®  Ua  .  sinacosg  (as)  q  (b3) 
sin  (0  *!■-!  r  YE  'I 


(aelj  +  22L2.  (avy 


sin  6  sin  a 
sin  (9  -ct  ) 


Jl  +  sinacosa.  (ds)j  „  0  (Bh) 


The  subscript  I  denotes  the  left-running  characteristic  and  the  subscript 
J  the  right-running  characteristic.  It  is  usually  more  convenient  to  work  in 
terms  of  Mach  Number,  M,  and  flow  angle,  9,  and  the  above  equations  can  be  transformed 
to  these  two  variables  through  the  following  relations: 


sin  a  = 


and  for  adiabatic  flow 


Using  these  relations,  and  the  corresponding  expressions  for  cosCX  and 
tan  a,  equations  (Bl)  through  (b4)  take  the  following  form: 


sin  9  V  m  -  1 
cos  -  1 


■f  cos  6 
-  sin  6 


sin  eVtf  -  1  -  cos  9 

cos  oVm  ■■  1  +  sin  9 


(a0)i  + 

V  t^-i  / 

YBlf 

ilFJ 

(d0)T  +  —  yV-O- 


sin  9 

0"V  bP-1  +  cos  9 


(d3)I  =  0 


_ sin  6 _ 

sin9  -  cos  0 


(BIO) 


♦  <«•>,  -® 

rai  1 
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The  numerical  solution  of  these  equations  covers  the  entire  flow  field 
through  a  combination  of  three  bacic  solutions  -  the  interior  point,  boundary'’ 
point,  and  shock  point. 

1.  Interior  Point  Solution 


Given  M,  0,  and  s,  as  well  as  the  coordinates  x  and  y,  at  two  initial 
points,  the  solution  can  be  extended  to  a  thirl  point  located  at  the  intersection 
of  the  characteristics  of  opposite  families  from  the  original  two  points.  Actually 
the  characteristics  are  curved  for  a  non-uniform  flow  .field,  but  are  approximated 
by  straight  lines  between  successive  points  with  a  slope  based  upon  the  average 
of  the  flow  properties  at  each  point .  Thus ,  for  accuracy,  the  size  of  the 
characteristic  net  should  be  kept  relatively  small. 


For  the  step-by-step  solution  to  this  "marching"  problem,  equations  (B7) 
through  (BIO)  must  first  be  put  into  finite  difference  form.  Using  the  notation 
in  the  above  figure. 


(Bll) 
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X3  X1  cos  Qj  “1  +  si-n  ®j 


(e  m'VS'1  1  (M3  -  V  t 
3  8  +  % 


(s3  -  s2)  ’  0 


778\/Mt  ~  1 


YK^ 


(e3  -  91>  + 


V1^  ~  1 
(i  +  i§) 


sin  0T 


_  .C5 -  _  ~Cr~+r^ 

sin  9j\/^  ”  1  +  cos  Qj  J  J 


(Mg  -  Mj_) 


r2^ 


sin  9n 


sin  ejVI'V  ‘  1  "  cos  °J 


(B13) 


2(r3-ri) 


tr3+V 


778-^17 

TKif 


(s3  -  sl)  "  0 


(BlH) 


and  Mj.  denote  average  values  along  the  I  characteristic  between  points  2 
and  3?  whila  Qj  aud  Mj  are  average  values  along  the  J  characteristic  between  1  and 

3. 

An  additional  equation  is  required  to  permit  a  solution  for  the  five  unknowns 

Xg,  r^j  and  9^  and  s^.  By  appropriately  interpolating  between  points  1  and  2  to 

locate  the  streauline  passing  through  point  3?  s_  can  be  expressed  in  terms  of  s.. 

49  3  1 

and  Sg.  Following  the  approach  of  Shapiro  ,  and  again  referring  to  the  last  sketch, 
normals  to  the  streamline  through  3  can  be  drawn  from  points  along  the  characteristics 
Since  the  entropy  can  change  only  along  the  normals,  the  change  in  entropy  along  the 
I  characteristic  can  be  expressed  as 


(as)j=  ^  (an^ 


Page  4 


♦ 

4 


where  (dn). 


I 

sin 
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sin  aT  , 

1  (dr),  -  — - -1 

da  9^*? 

can  be 


(e,  +  a. 


x  y  j 

The  normal  entropy  gradient  can  be  approximated 
of  points  1  and  2  divided  by  the  distance  between  the  1 
measured  along  the  normal  to  the  stream  line  through  3; 


~v 

+  cos 


sin  9£\J  1^-1  x 

approximated  by  the  difference  in  entropy 
between  the  streamlines  through  1  and  2 
Lne  through  3-,  or 


ds  ^ 
dn 


S1  -  s2 

%  +n2 


(x,  - 


x^)  sin  Hj  (x^  -  Xg)  sin 


(B16) 


cos  (9 


(® T  -  ar) 


+  a-j.) 


ables 


cos  (9^ 

Therefore,  substituting  equation  (Bl6)  into  (B15),  transforming  to  the  vari- 
1  and  9,  and  writing  the  result  in  finite  difference  form  gives: 

r  si  -  s? 


s3  -S2 


-  X  ) 


cos  9j 


X  _  ^  _  D  mm  ^ 

1  +  sin  9j  cos  M^  -  1  -  sin  9 ^ 


_ r3  ~  V2 

sin  -  1 

Equation  (B17)  added  to  equations  (Bll)  through  (Bl4)  constitute  the 
set  of  simultaneous  equations  to  be  solved  for  x^,  r^,  M^,  Q^,  and  s^« 

This  solution  proceeds  in  the  following  manner: 

a.  Setting  Mj.  =  Mg,  Gj  =  9g,  M.  =  M^  9..  =  9.^  equations  (Bll)  and  (B12) 
are  solved  simultaneously  for  and  r^. 

b.  Equation  (B17)  is  solved  for  s^. 


cos  9^. 


(B17 


c.  Equations  (B13)  and  (Bl4)  are  solved  simultaneously  for  Ma  and  9 


MU  +  It,  9„  +  9„  1^  +  M3  9,  +  9 


Now  setting  Mj.  =  -2  g — — ,  9^.  =  — — -,  Mj  = 


1  3 

and  9r  = — ~, 
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Equations  (Bll)  and  (B12)  are  again  solved  for  new  values  of  x^  0n^ 

Steps  b,  c,  and  d  are  repeated  until  successive  values  of  satisfy  the 
following  criterion: 

KM3>n  -  I  <  -0001 

The  total  pressure  at  point  3  is  a  function  of  the  entropy  at  3  and  the 
total  pressure  and  entropy  at  some  reference  station 


P,  =  P,  Exp 
*3  *o 


778(So  -  s3)  -] 

T> 


(BIB) 


For  the  LITVC  flow  field  solution,  P^  is  taken  as  the  chamber  pressure 
and  sq  is  set  equal  to  zero.  The  static  pressure  can  now  be  calculated  from  the 
isentropic  pressure  ratio  based  upon  P^,  which  is  constant  along  the  streamline 
through  point  3. 

Y  ,  2  _r_ 

P3  -  Pt3  (1  +  M3  )  1_Y  < 


2  3 


2.  Boundary  Point  Solution 


(B19) 
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In  the  L-ITVC  flow  field  solution,  the  pressure  must  he  calculated  along 
the  solid  boundary  formed  by  the  vapor  body  and  forward  separation  region. 

» 

Knowing  the  flow  properties  at  two  points,  along  an  upstream  I  characteristic, 
t  and  x,  r,  9,  and  s  at  the  solid  boundary  point,  the  Mach  number  at  1  can  be  calcu¬ 

lated  by  the  following  procedure. 

M  +  M, 

a.  The  initial  values  for  and  are  set  equal  to - - -  . 

_  M  +  m4  _  02  +  \  _  m  +  M  _  e  +  0 

b.  Setting  Mj  =  . . . -,  0J  =  •  —  g . -,  Mj  =  -  ±-  —^,  and  0j  = . ■  -, 

''  point  3  is  located  by  solving  equations  (Bll)  and  (B12)  as  in  the  interior  point 

solution . 


c.  M^,  ©g  are  calculated  by  linearly  interpolating  along  the  I 

characteristic  between  points  2  and  4,  i.e. 

(x  -  x_) 

VV  J  (m4-“2> 

d.  Equation  (Bl4)  can  now  be  solved  directly  for  M^,  and  steps  b,  c, 
and  d  repeated  until  successive  values  of  satisfy  the  criterion 

I  (“I'n  -  (Vn-l  I  <  -O001 

The  static  pressure  can  then  be  calculated  with  equations  (Bl8)  and  (B19) • 
3.  Shock  Point  Solution 
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For  a  given  Mach  number  before  an  oblique  shock,  M^,  and- a  given  flow 
deflection, 6 ,  through  the  shock,  the  Mach  number  behind  the  shock,  M^,  the  shock 
angle,  e ,  with  respect  to  M^,  and  the  entropy  rise  across  the  shock,  -  s^,  can 
be  calculated  with  the  following  relations: 


cot  5  =  — g—g - “  1)  tan  e 

Mx  sin  £  -  1 


1 _  f  Y  +  1  tan  e 

£  -5  )  L  2  tan  (& 


sin( £ -6 


XzJ,  ‘ 


(B20) 


(B21) 


sa  -  \  ■  y^i  10  f  ffi  Mis  si"2  6  Tri  ] 


(B22) 


+  P-  In  [ 

Y-  1  L 


r- 1 ' 
+ 1 . 


(Y  +  l)!^2  sin2e  Y 


J  CHARACTERISTIC  iN  FRONT  OF  SHClXP 


oHtv  n  s  N 


9Ai  A 


^  >v  '' 


'I  CHARACTERISTIC  INFKWTf'W 


'I  CHARACTERISTIC  BEHIND  SHOCK 


J  CHARACTERISTIC  Pt-IKO  ->hC« 
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To  calculate  the  shock,  point  B,  the  flow, properties  must  be  known  at  point  E 
in  front  of  the  shock,  at  point  A  both  in  from  of  (Al)  and  behind  (A?.)  the  shock'/'' 
and  at  point  D  on  the  same  J  characteristic  behind  the  shock  as  A2-'  The  following 
procedure  is  then  used: 

a.  The  imaginary  point  F  is  calculated  with  the  interior  point  solution 
using  the  properties  at  E  and  Al- 


f* 


b.  For  6  =  ©A2  "  ®ai  =  ^Al’  shock  angle  at  A,  E^, 

with  equation  (B20)  using  an  iterative  procedure. 


is  calculated 


c.  Point  B  is  located  at  the  intersection  of  the  J  characteristic  between 
E  and  F  and  the  shock  extended  from  A  at  the  angle  (eA  +  0^)  with  the  horizontal. 

d.  At  B,  the  flow  properties  in  front  of  the  shock  (Bl)  are  calculated  by 
interpolating  between  E  and  F. 


e.  From  Equation  (B20)  with  = 

©B2  =  +  ^  b  •  is  catculated  from  equation  (B2l)  and  from  equation  (B22) . 

~  —  ^Pi2  +  —  ®A2  + 

f.  For  Mj.  =  Mgg,  ©j  =  ©B2,  Mj  =  - 2 - ,  and  Gj  = - 2 - >  equations 

(Bll)  and  (B12)  are  used  to  find  point  C. 

«* 

g.  The  flow  properties  at  C  are  found  by  interpolating  between  A2  and  D. 

equations  (Bll)  and  (B12)  are  again 


I.  (“c’n  '  <MC>n-l  I  <  •°001 


h.  For 


.  i4B2+MC  *  9B2t0C 


*1 


QI  = 


used  to  find  a  new  point  C. 

i.  Steps  g  and  h  are  repeated  until 


and  £  =  E^,  6  g  is  calculated,  and 


i 

t 


Page  9 


TM-16-SR0,  Appendix  B 


3 .  Equation  (B13)  is  used  to  calculate  a  new  value  for  0__,  which  is 

Bd 

averaged  with  the  last  value  to  calculate  a  new  6  =  ©_0  -  ©__ . 

D  u  Bd  BJ. 

k.  The  shock  angle  at  B,  £3  ,  is  calculated  with  Equation  (B20)  for 

"l  =  liBl  ?  -  V 

l.  A  new  point  B  is  located  by  extending  the  shock  from  A  now  at  the 

£A  *r  £B\ 

angle  ( - - )  +  9^  to  the  horizontal  and  intersecting  the  J  characteristic 

between  E  and  F.  ' 

m.  Steps  d  through  1  are  repeated  until  successive  values  of  satisfy 
the  criteria  |  (Mgg)n  -  (Mgg^n-l  I  <  .0001. 

Once  the  flow  properties  at  B2  are  found,  they  are  used,  along  with  the 
properties  at  D,  to  find  a  new  interior  point  and  continue  the  solution. 

The  three  basic  point  solutions  just  described  have  been  combined  into  a 
computer  program  for  the  IBM  709^,  which  takes  less  than  1  minute  per  case  for  a 
typical  flow  field.  The  input  required  for  this  program  is  M,  0,  s.  x,  and  r  at 
a  number  of  points  in  front  of  the  shock,  and  0,  s,  x  and  r  at  points  along  the 
lower  boundary.  The  lojie  followed  by  the  program  is: 

•-  -  a.  Calculate  the  shock  point; 

b.  Continue  the  solution  by  calculating  interior  points  along  the  J 
characteristic  from  the  shock  point  until  the  lower  boundary  or  last  I  character¬ 
istic  is  reached; 

c.  Calculate  the  next  lower  boundary  point,  if  there  i3  one,  and  the 
interior  point  at  the  intersection  of  the  I  characteristic  from  the  boundary 
point  with  the  J  characteristic  being  computed. 
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d.  Repeat  a,  b,  and  c  until  all  the  initial  points  in  front  of  the 


shock  are  used. 


Figure  B1  compares  the  shock  shape  and  pressure  distribution  calculated 
by  this  program  with  experimental  data.  Very  good  correlation  has  been'obtained. 
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(Ref:  NASA  TN  D-1436,  pp  28  & (30)  '  \ 
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- -  Taken  from  Schlieren  Photographs  ' 

»  Calculated  shock  point 
■  *  ''  Calculated  separation  shock  point 

assuming  separated  region  to  act 
as  a  solid  boundary  ■* 


1  M*  =  4.98 


Flow  Separation  Region- 
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Comparison  of  Calculated  Shock  Shapes  and  Pressure  Distributions 

with  Experimental  Data 


*  Figure  1 


